Spin-electric stripes: Electric voltage induced by spin currents 
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At each of the boundaries of the two-dimensional (2D) rectangular conductor parallel to the 
electric current there arises a stripe with an electric field transverse to the current and a 100% 
electron spin polarization. The two stripes have opposite spin orientations and opposite directions 
of electric fields. The magnitudes of the fields, directly related to the spin current if the spin 
relaxation is negligible, are the same. The periphery stripes are separated by a center-stripe, in 
which the magnitude and direction of the electric field depend on the ratio of the skew scattering 
and side jump spin currents. The spin polarization is zero on the centerline and reaches ±1 at the 
boundaries between the central and periphery stripes. Weak relaxation of the z— component of spin 
normal to the 2D plane modifies the magnitudes of the spin polarization and fields, with ±1 spin 
polarization persisting at the edges of the sample. Favorable experimental settings, in which the 
electron spin relaxation of the z— component of spin is suppressed but the spin current is not, are 
discussed. 

PACS numbers: 72.25Dc, 72.25Rb, 71.70EJ, 85.75-d 
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Introduction When an electric current is flowing in a 
conductor, the spin-orbit interactions result in the asym- 
metric scatteringiJ, an accumulation of spin polarization 
at the sample boundary^], or an average electron spin 
polarizatiorPEl xhe accumulation of spin polarization 
at the boundaries in the presence of electric current has 
been the subject of recent experimentaStiS aimed at the 
observation of the spin Hall effects discussed theoreti- 
cally over the year^^UiHlSl xhe conventional wisdom is 
that no ampermeter or voltmeter exists that can directly 
measure the spin Hall current or spin Hall voltage (the 
voltage and current signals due to charge carriers of op- 
posite spins cancel each other). Therefore these exper- 
iments address spatially-dependent electron spin polar- 
ization S(r). S(r) is generally defined by several effects 
besides the spin Hall current. The total spin density 
is determined by the spin polarization of electrons in- 
duced by the electric current, by the spin relaxation and 
spin diffusion, and by the spin precession in the effec- 
tive momentum-dependent fields (for example, the fields 
due to the Dresselhaus or Rashba coupling) and the ex- 
ternal magnetic field. In general, there is no continuity 
equation for the spin current in contrast to the continuity 
equation for the charge current. Measurement of S(r) (or 
time-dependent S(r,t)) is not capable of distinguishing 
between the several contributions to the net spin density. 
Furthermore, the inhomogeneous spin density arises due 
to the boundary effecta^SEl]. This happens even when 
the bulk spin Hall conductivity completely vanishes, e.g., 
in systems with o nly in trinsic spin-orbit interactions lin- 
ear in momenturrP32^i25l_ Moreover, the charge current 
passing through the sample carries its own magnetic field, 
whose directions at two edges of the conducting film or 
the 2D plane are opposite to each other. The polarization 
of electrons due to the Zeemann interaction caused by 
this magnetic field leads to generation of opposite spins 
at two edgea^ESl. In semiconductor materials with a large 
g-factor, this contribution into spin polarization can ex- 




FIG. 1: Schematic picture of electron density distribution in 
a 2D sample of width 2a, with two periphery stripes with op- 
posite electric fields transverse to a fiowing current separated 
by a zero-field central stripe with width 26. The picture cor- 
responds to case I of Sec. lA. In the periphery stripes, only 
electrons with one projection of spin (-|- or -) contribute to the 
electric fields. In the center-stripe the gradients of density ex- 
ist for both spin projections, and the transverse currents due 
to an applied external field for each of the projections are 
cancelled by the corresponding diffusive currents. 

ceed the contribution to the spin polarization caused by 
the spin current. Such spatially dependent polarization 
is not related to the spin Hall current but is indistin- 
guishable from it experimentally when the spin density 
is studied. Also, the average electron spin polarization 
in the presence of electric current observed in^ '^^, while 
contributing to the overall spin density, is not related 
to the spin Hall effect in any direct way. The question 
arises: is it possible at all to find an unambiguous sig- 
nature of the spin Hall current showing no or negligible 
contribution due to the other effects? 

In this paper, we predict that in two-dimensional (2D) 
systems with spin-orbit interactions that conserve an 
electron spin projection perpendicular to the 2D plane. 



the electric current results in two stripes with uniform 
electric fields of opposite signs and equal magnitudes at 
two opposite boundaries of the sample. The electric cur- 
rent flows along the stripes, and the electric fields are 
perpendicular to the flowing current, and are directly re- 
lated to the spin Hall current. In the stripes, electrons 
are fully spin polarized, with opposite spin orientations 
at opposite sides of the sample. The two stripes are sep- 
arated by a third stripe in the central region of the sam- 
ple. On the center-line, the electron spin polarization is 
zero. The picture of the spin polarization in the cen- 
tral stripe is almost like in the conventional spin accu- 
mulation, except that when the spin polarization grows 
from the center-line to the boundaries between the cen- 
tral and the periphery stripes (or falls in the case of the 
opposite spin direction), it reaches maximal ±1 value at 
that boundary, and cannot grow (fall) further into the 
periphery stripes. The magnitude and direction of the 
electric field in the central stripe depends on the relation 
between microscopic mechanisms of the bulk spin Hall 
current. For exampl e, if t he only spin current is due to 
side-jump-like effectP^ESl, then the electric field in the 
central stripe vanishes, as in Fig. [T] In contrast, the 
skew scattering spin current results in the electric field 
in a central stripe, which is linear in transverse coordi- 
nate, vanishes at the center-line and is quadratic in the 
current flowing through the sample. The symmetry with 
respect to the center line characterizes the system in the 
absence of an external magnetic field only. In the pres- 
ence of the external magnetic field, the stripes become 
asymmetric. We will discuss various experimental situa- 
tions and their influence on the transverse electric field 
profile in stripes. 

The origin of the different roles of side-jump-like and 
skew scattering spin currents in the stripes is their differ- 
ent dependence on the density of contributing electrons 
and on the spin density and polarization. Side-jump-like 
spin currents are independent of the spin density and de- 
pend only on the total density, in contrast to the skew 
scattering spin current, which depends on both densities. 
If the spin polarization varies in space, the spin current 
can change its magnitude and even its direction. This 
may result in a non-trivial spatial disribution of electric 
field in stripes. 

The transverse electric fields result in a potential differ- 
ence between the central stripe and the edges of the sam- 
ple. The example of a potential profile that corresponds 
to Fig. [T] is shown in Fig. [2J Measuring the voltage in an 
experiment similar to^ is a straightforward method of 
observation of the spin-electric state. Positioning probes 
to measure voltage caused by the electric fields in the 
periphery stripes would constitute a direct electric mea- 
surement of the spin currents. Furthermore, if a single 
periphery stripe becomes a part of a separate electric cir- 
cuit with the current flowing, it is possible to transfer its 
giant spin polarization for the purpose of applications. 

Spin-electric stripes are accompanied by the charge re- 
distribution in the sample. The electric fleld in the pe- 
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FIG. 2: The Hall direction potential distribution across the 
sample (case I Sec. lA). Slopes of potential curves on the 
right and on the left correspond to opposite electric fields in 
two periphery stripes, and are directly related to the spin 
Hall effect and are linear in the flowing electric current. The 
underlying parameters are discussed in Sec. IV. 
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FIG. 3: Cartoon showing the physical origin of the periphery 
stripes with opposite electric flelds due to the spin current in 
the absence of z-projection spin relaxation. When spin polar- 
ization reaches ±1 at the boundaries of the periphery stripes, 
this spin polarization makes the spin current equal to the 
Anomalous Hall current, which leads to charge accumulation 
near boundaries of the stripes. The Coulomb interactions pre- 
vent electrons from mounting walls of the sample and results 
in electric fields. 

riphery stripes is linear in the externally applied electric 
fleld, and the flowing current. These are "true" elec- 
tric flelds an d cha rge distribution, which do not require 
nonlinearitiea^^'^fil or pre-set, artiflcially created macro- 
scopic inhomogeneity of charge density in the sample^. 
Charge carriers with both spin projections contribute to 
the macroscopic electric field and charge density distri- 
butions. 

Furthermore, we find that spin-electric stripes arise 
even in the presence of spin relaxation, if the latter is 
sufficiently weak. As we demonstrate, macroscopic con- 
sequences of the spin current differ in systems with strong 
and weak spin relaxation. In systems with strong spin re- 
laxation, it is the latter that balances the spin current, 
keeps the total concentration of charge carriers constant 
across the sample, and regulates spin accumulation at 
the boundaries. The spin current induced by an external 



electric field in the presence of spin relaxation results in a 
spatial dependence of the spin density. Correspondinly, 
a diffusive spin current arises, which balances the ex- 
ternally induced spin current, particularly canceling the 
total spin current at the surface. At strong spin relax- 
ation, stripes do not arise. Once spin relaxation is not 
present, or is weak (the criteria for the spin relaxation 
strength will be discussed below), it is no longer capa- 
ble of balancing the spin current. Spin polarization and 
density are limited only by their maximal (minimal) val- 
ues, which are reached in the periphery stripes. Here 
the spin polarization is a sourse of the Anomalous Hall 
effect (AHE), which results in a charge build up. This 
buildup is opposed by the Coulomb interactions and elec- 
tric fields, somewhat similar to the usual Hall effect, but, 
of course, with opposite electric fields in two stripes, in 
contrast to the usual Hall effect. For example in the case 
shown in Fig. [2] that corresponds to side-jump-like cur- 
rents only contributing to the spin Hall effect, the poten- 
tial profile caused by transverse electric fields is aimed at 
keeping electrons in the sample. Were there no Coulomb 
interactions, charges with opposite spins as a result of 
multiple scattering and spin current would stream to- 
wards the opposite boundaries and mount the walls of 
the sample, leaving the middle of it free of charge. The 
Coulomb interaction prevents this (Fig. ^ and results in 
the electric field and spin-electric stripes. We find that 
because the sample breaks in the periphery stripes and 
a central stripe, and thus becomes non-uniform, other 
curious spatial distributions of the electric field and po- 
tential in the stripes are possible, depending on the ratio 
of the side-jump-like and skew scattering spin currents, 
and their variation across the sample. 

The ideal picture in the absence of relaxation of the 
z— component of spin is not much affected by weak spin 
relaxation. We show that in this case, spin polarization 
and the electric field in the periphery stripes experience 
a small change from the periphery towards the center 
region. 100% spin polarization is still present at the 
edges of the sample. The division of the sample into 
three stripes remains meaningful in the presence of weak 
spin relaxation: at the boundaries between the periphery 
stripes and the central stripe, the electric field is a steep 
function of the transverse coordinate. 

The state of the matter that we describe is different 
from that discussed by Dyakonov and PereP or Hirsclr^. 
The distinction is both on macroscopic and microscopic 
levels. On the macroscopic level, no voltage due to the 
spin Hall effect between the center of the sample and its 
sides arises in these studies (the voltage probe suggested 
by Hirsch is in the longitudinal direction, is quadratic in 
the anomalous Hall resistance and constitutes a signal 
that is three orders of magnitude smaller than the trans- 
verse voltage we discuss). Furthermore, the phenomenon 
we predict here arises due to the spin Hall current in 
the presence of significant spatial separation of electrons 
with opposite spin projections. Moreover, discussing the 
spin-electric state in the presence of spin relaxation, we 



will see where the traditional spin accumulation physical 
picture fails and the spin-electric stripes emerge. On the 
microscopic level, the state we suggest arises due to elim- 
ination or suppression of the part of the spin-orbit inter- 
action responsible for the relaxation of the z— projection 
of spin, while the spin-orbit interactions that produce the 
spin current are kept intact. This physics is described by 
the Hamiltonian different from those studied in^ ^■^. 

We consider macroscopic properties of the spin Hall 
system in terms of an ideal picture of the 2D metal. 
We also discuss the role of the effects of inhomogeneous 
distribution of impurity centers in semiconductor het- 
erostructures, and show that they do not alter the prin- 
cipal macroscopic observables in any essential way. 

Spin-orbit interactions that conserve the electron spin 
projection transverse to the 2D plane arise for generic 
spin-orbit interactional 
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H = 



2m 



[/(r) + aa^ ■ [k x VrC/(r)] 



(1) 



where k and r are the two-dimensional wavevector and 
coordinate, m is the effective mass, tJz is the z— Pauli ma- 
trix, C/(r) includes potentials due to impurities and due 
to an external electric field, and a is the spin-orbit con- 
stant (see Appendices B and C). The remarkable prop- 
erties of this Hamiltonian in 2D metals, containing just 
one Pauli matrix, because of the presence of only in-plane 
electron moment a, hav e been noticed in studies of weak 
localization effect^^SMOl Spin accumulation in this system 
has been also calculatecP^. A similar situation with con- 
served spin can arise in materials with intrinsic spin-orbit 
interactions of the Rashba and Dresselhaus type, but only 
for special crystallographic orientations or for equal mag- 
nitudes of the Rashba and Dresselhaus constanti^2M6|_ 
For spin-orbit interactions described by Eq.pl), no spin 
precession, no generation of average spin by the electric 
field and no spin relaxation of the z-component of spin 
takes place. Therefore, the spin current Js is unam- 
bigously defined. 

Even the simplest spin-orbit interactions Eq. M lead to 
several microscopic contributions to the anomalous Hall 
effect and the spin Hall curren t. In this work, we use 
a universally accepted approactPSBIHSl] for the spin cur- 
rent, and calculate it in cases favorable for observation of 
the spin-electric stripes. The role of the remote doping 
by donors and doping in the quantum well by donors and 
acceptors is studied. By using various degrees of compen- 
sation of donors by acceptors in the quantum well, it is 
possible to eliminate skew scattering altogether. In this 
respect 2D structures with remote doping are a more fa- 
vorable experimental setting compared to 3D samples^. 

Simple 2D metal systems, in which the electron spin- 
orbit scattering is described by the Hamiltonian (fTl) , have 
not been acheived technologically. We therefore turn our 
attention to semiconductor heterostrutures. In order to 
find a feasible experimental setting, we embark on a com- 
prehensive study of the spin-orbit interactions in quan- 
tum wells. We shall demonstrate that there can exist a 



class of material systems, in which spin-orbit interaction 
leading to the spin current Js exists, but spin-orbit in- 
teractions resulting in spin relaxation of the z— projection 
of spin vanish. In these systems, the intrinsic spin-orbit 
coupling of the Dresselhaus type and strain-induced in- 
trinsic spin-orbit interactions are removed by using a spe- 
cial crystallographic orientation of the 2D plane. Intrin- 
sic spin-orbit interactions of the Rashba type are removed 
using symmetric potential confinement of the quantum 
well, leaving only the residual Rashba interactions asso- 
ciated with large scale fluctuations of density of impuri- 
ties located in symmetrically positioned (with respect to 
the quantum well) ^-doping layers^. In our considera- 
tion of the spin-orbit effects for 2D electrons, 3D short- 
and long-range random potentials are taken into account. 
The role of the boundary conditions and the role of the 
difference of parameters in the barriers and the quantum 
well is investigated. We demonstrate that by selecting 
heterostructure materials, it is possible to eliminate this 
residual Rashba interaction, as well as spin-flip scattering 
due to impurities of the quantum well. The only source of 
spin relaxation in these systems is then extremely small 
interactions of the free electron spins with the spins of nu- 
clei. It turns out that not only are such materials theoret- 
ically possible but also there already exists an electronic 
system in InAlAs/InP/InAlAs double heterostructure, in 
which it looks feasible to observe the spin-electric stripes 
caused by the spin Hall effect. 

The manuscript is organized as follows: In Sec. I, 
we present the discussion of the stripe structure due to 
spin Hall current in the 2D metal. Sec. lA deals with 
a macroscopic equations and properties of the state in 
the absence of spin relaxation, IB describes the charge 
distribution corresponding to the stripes, IC describes 
spin-electric stripes in the presence of magnetic dield per- 
pendicular to the 2D plane, I.D discusses a macroscopic 
picture in the presence of non-uniform doping, and Sec. 
IE discusses microscopic mechanisms of the spin current. 
In Sec. II we discuss the spin-electric state the presence 
of spin relaxation. Sec. HA gives macroscopic equations 
in the presence of spin relaxation, IIB shows how the 
traditional spin accumulation picture stops working and 
a spin-electric state emerges, IIC gives solutions of the 
macroscopic equations for the spin polarization and elec- 
tric fields and presents a criterion of weak relaxation of 
the z— component of spin. Sec. HI discusses a strength 
of spin-orbit interactions and spin relaxation in various 
electronic systems and feasibility of spin-electric stripes. 
In particular, IIIB discusses spin-orbit interactions in 
the quantum wells, and how they can be eliminated or 
strongly suppressed in order to observe the spin-electric 
state, IIIB. 3 describes impurity spin-orbit interactions in 
2D systems, IIIB. 4 discusses the residual Rashba intrin- 
sic spin-orbit interactions due to large-scale fluctuations 
of doping potential in heterostructures, IIIB. 6 is devoted 
to the Dyakonov-Perel spin relaxation for non-uniform 
spin-orbit interactions, which is the case for the residual 
Rashba interaction, IIIB. 7 analyses the role of intrinsic 



spin-orbit coupling for the non-equilibrum and equilib- 
rium spin currents in [110] grown quantum well. Sec. IV, 
which discusses feasible experimental setting for observa- 
tion of the spin-electric state in InAlAs/InP/InAlAs dou- 
ble heterostructure, is followed by a conclusion. There 
are three appendices in the paper: Appendix A discusses 
details of the solution of the 2D electrostatic problem 
that defines self-consistently the distribution of charge 
density and electric fields in the sample. Appendix B 
discusses spin-orbit constants in various systems, includ- 
ing the role of boundary conditions in heterostructures, 
and treats spin-orbit constants in quantum wells in the 
III-V materials. Appendix C discusses the effect of dif- 
ferent materials parameters in the quantum well and the 
barriers on spin-orbit constants responsible for relaxation 
of the z-component of electron spin. 



I. SPIN-ELECTRIC STRIPES IN THE ABSENCE OF 
SPIN RELAXATION 

The principal goal of this section is to solve macro- 
scopic constituitive equations, a macroscopic electro- 
static problem in the 2D electron system, and to find 
the electric field and the electron density profile in the 
presence of the spin current. We will also discuss the mi- 
croscopic theory of spin current aiming at experimentally 
relevant situations, and describe the stripe state in the 
presence of an external magnetic field perpendicular to 
the plane of the 2D gas. 



A. Constitutive equations of the system 

The system is assumed to be a rectangle in the xy 
plane, as shown in Fig. [4j with a width 2a in z-direction 
and the length of the sample between the electrodes in 
j/-direction L 3> a, considered infinite. The boundary 
conditions for this configuration is that there be no cur- 
rent flow in x-direction out of the sides of the device, and 
no electric field parallel to the contacts in j/-direction. 
The electric fleld E = — V(/>, where is the electric po- 
tential, and the current density is j(n), where n = erie 
is the charge density and Ue is the electron density, are 
assumed to be independent of the y-direction, i.e., we are 
dealing with the strictly a 2D electrostatic problem. 

The fundamental constitutive equations for the system 
are 



j± = cr(n±)E + D{n±)Vrn± + 
crj_(n±)(E X z) + Dj_{n±){Vrn± x z). 



(2) 



where signs ± correspond to the spin z-projection, which 
is conserved for systems described by the Hamiltonian 
Eq. (IT]) , n± is the concentration of carriers with the corre- 
sponding spin projection, a(n±) is the longitudinal con- 
ductivity due to carriers with ± spin projection, cr±{n±) 
is the transverse conductivity, D(n±) is the longitudinal 




FIG. 4: Schematic picture of the spin Hall sample. 2D elec- 
tron is gas confined to a quantum well of width d. 

diffusion constant describing the current directed along 
the gradient of the carrier concentration, and D±{n±) 
describes diffusion in the direction transverse to the con- 
centration gradient. The longitudinal component of these 
equations expresses Ohm's law and corrections to the cur- 
rent due to spin polarization. The transverse component 
vanishes, because of vanishing currents at the boundary 
of the 2D plane and their continuity in the absence of 
spin-relaxation of the z— component of spin. Expressions 
for all the coefficients in these equations are obtained 
from the microscopic description. The Einstein relations 
hold for both longitudinal and transverse transport, and 
give relations between a{n±) and D{n±), and between 
cr_L(n±) and Di_{n±). 

Eqs.([2]) for the current of charge carriers with spin pro- 
jections "-I-" ("— ") contain carrier density, conductivity, 
and longitudinal and transverse diffusion constants cor- 
responding to the same sign of the spin projectiorP^]. The 
electric field E(x, y), which enters the equations for both 
spin projections, is defined, apart from the external elec- 
tric field, by all charges in the system, including carri- 
ers of both spin projections and impurities supplying the 
charge carriers (for doped semiconductors) or ions (for 
metals). 

We now describe the physics and find the solutions of 
these equations. As a starting point, we take a uniform 
sample. Assuming a certain sign o f a (the sign of a can 
differ from system to systenPSMl) ^ charge carriers with 
spin projection "-I-" leave the boundary 2: = a, and car- 
riers with spin "— " leave the boundary x — —a because 
of the spin current (Figs, [l] and p|. Furthermore, charge 
carriers with any transverse coordinate and spin projec- 
tion "+" attempt to move in the negative direction of 
x, towards —a, and carriers with the spin projection "— " 
attempt to move in the positive direction of a;, towards a. 
We now consider an ideal case with positive charges (due 
to impurities or ions) uniformly distributed over the sam- 
ple, in the conditions of the overall electrical neutrality. 
Then the sum n+ + n_ , which is the total charge carrier 
density n, in the first approximation can be considered 



constant with transverse coordinate x. 



n(x) — Ti^{x) + n_{x) — Const 



(3) 



because in such situation the Coulomb interactions gener- 
ally favor homogeneous distribution of carriers. However, 
the densities of electrons with each component of spin 
n^{x) and n^{x) are not uniform, and therefore coeffi- 
cients in Eq.S are coordinate-dependent. Furthermore, 
the spin polarization 



Sz{x) = 7T._|_(a;) — n^{x) 



(4) 



is not uniform. At the periphery of the sample a non-zero 
spin polarization arises defined by the sign of the spin- 
orbit constant defining the spin current. In conditions 
of a uniform total density, the spin polarization at the 
centerline (centerplane) of the sample, on the symmetry 
grounds. 



5.(0) =0. 



(5) 



We note that by our definition, Sz{x) has the same 
dimensionality as the charge density n. We will primarily 
use Sz{x) and call it, for simplicity, the spin density. The 
polarization is given hy p = Sz{x)/n, and electron spin 
polarization is Tip/2. The spin density in the strict sense 
is h.Sz{x)/2e, where e is the value of the electron charge. 

In order to solve the constituitive equations ([2|, we 
need the dependences of the coefficients in these equa- 
tions on density of charge carriers (with a given projec- 
tion of spin). We assume that the conductivity of the 2D 
gas is defined by its Drude value 



c("'±) — eTn±/m = ^n±, 



(6) 



where r is the transport relaxation time, and /i is the 
electron mobility. We note here that for two indepen- 
dent electron species with up and down projections of 
spin, spin polarization effectively leads to two quasi- 
Fermi energies^. In terms of the degree of electron po- 
larization p, these energies in two dimensions are given 

by 



e| = 6^(l±p), 



(7) 



where e^p is the Fermi energy for unpolarized electrons 
(p = 0). The charge densities of electron species with up 
and down spin projections are 



n± 



em 

2TTh 






(8) 



In a degenerate electron gas, generation of electron-hole 
pairs (an electron above quasi-Fermi-energy, and a hole 
below it) for minority species takes place at a different 
energy compared to majority electrons. This happens 
because the Pauli principle allows for unoccupied states 
for minority species in the vicinity of the lower quasi- 
Fermi level. Correspondingly, the Drude conductivity 
is a sum of conductivities of up and down species, with 



each partial low-temperature conductivity defined by the 
quantities taken at the corresponding quasi-Fermi level. 
Other kinetic coefRcients defined by the quantities at the 
Fermi level for unpolarized electrons, for two independent 
electron species are defined similar to conductivity. The 
longitudinal diffusion coefficient for each of the electron 
species is given by 



the spin current and in Eq.(12) for the charge current 



D{n±) — ttTi n±T/em = C"-±- 



(9) 



Here we assume that r is energy independent, and has the 
same value for both spin projections. We also note that 
the expression for the longitudinal diffusion coefficient 
explicitly assumes the 2D case. The coefficients ( and /i 
are independent of densities. The generalizations taking 
into account the spin projection dependence of r, ^ and 
/i are straightforward. 

We now observe that in the first approximation, there 
is no gradient of densities for up and down spin pro- 
jections in the longitudinal direction y. Such gradients 
appear only in the second order in the spin-orbit con- 
stant (one order comes from the transverse spin current, 
and the other is required to generate the coordinate de- 
pendence of density in the y-direction, when spin-orbit 
scattering transforms the transverse spin current into the 
longitudinal one). The effects of the term proportional 
to D±{n±) are cubic in the spin-orbit constants, because 
Dj_{n±) is connected by the Einstein relation to a± and 
is linear in the spin-orbit constant. Therefore, such ef- 
fect are small and from now on we neglect terms with 
-D±(n±). 

For the solution of Eqs.(l2]) wc also need the coordinate 
dependence of the transverse conductivity. As we shall 
see in Sec. IE, the transverse spin conductivity in rel- 
evant systems is defined by two principal contributions: 
one is linear in the electron density (this contribution 
originates from the effects of the side-jump group) and 
the other is quadratic in density (and comes from skew 
scattering) : 



cr±in±) = ±{riin± + m^^-) 



(10) 



To make observations on the physics of the phe- 
nomenon, it is instructive to write down the identities 
for transverse spin and charge currents by adding and 
subtracting Eqs.(l2| and taking into account that these 
currents vanish, J^ = because of the boundary condi- 
tions and continuity in the absence of spin relaxation: 

^iS^E^ + ^S^. + (2%" + min^ + Si)) ^ = 0(11) 
finE, -t- ||:5f + im + mn)S,Ey = 0, (12) 



are different in general. Only in the case of fully polar- 
ized electrons {S^ — i^T-), when the values of the spin 
current and the AHE current are the same, the combina- 
tions of skew scattering and side-jump-like terms in both 
equations become ±(771 -I- ri2n) ■ n and are the same. At 
small Sz, the spin current and the AHE current differ, 
and that will have important consequences. We under- 
score that the spin current depends on spin density via 
the skew scattering contribution only. In the case of vary- 
ing spin polarization, the value and even the direction of 
the spin current can vary across the sample. 

Vanishing transverse spin current Eq.(ll| is a result 



of the balance of the spin mobility current, which arises 
if S'2 7^ and E^. 7^ are generated due to spin-orbit 
interactions, a diffusive spin current, which arises if Sz 
is not uniform, and the bulk (sourse) spin Hall current. 
Vanishing transverse charge current ( 12 ) is a balance of 
mobility current steming from E^, diffusive current aris- 
ing due to non-uniform n^ + n'^_ = (S*^ -|- n^)/2, which 
is equivalent to non-uniform Sz at unifom n , and the 
AHE current. These identities are non-linear equations 
for Sz and E^, and solving them is not straightforward. 
Instead, to find a solution of the problem, we return to 
the constituitive equations ^ for j"^ —Q. 

We first solve these equations for the case when both 
n+(a;) and n-[x) are non-zero: 



n^{x) > 0; n_{x) 



>G. 



(14) 



We then divide Eq.([2| for j+ = by n+ and Eq.S for 
j^ = by ri_ , and obtain 

f^E^ + Ci^n+ + {jn + V2n+)Ey = (15) 

fiE^+(:^n_~{m+V2n-)Ey^Q. (16) 

Since we neglected terms with D±{n±), in these equa- 
tions Ey is an externally applied electric field that gen- 
erates a current in the longitudinal direction, and Ex{x) 
is the transverse electric field that we have to find along- 
side with the charge density n+ (n_ = n — n^). Adding 
Eqs.(15l and (16), and taking into account Eq.(l3]), we 
obtain the equation for the transverse electric field 



E^, 



'ri2Sz{x)Ey 



(17) 



Thus, the transverse electric field is defined by Sz and 



the skew scattering constant 772, ff Eqs.(14| are satisfied 



The skew scattering term in Eq. (17) arises as a specific 



where n\ - ?ii = Sz-n, and the strusture of the 772 -term for of spin polarization 



linear combination of the AHE and spin Hall currents. 
Subtracting Eq. ( 16 ) from Eq. ( 15 ) , we obtain the equation 



comes from (n_ 



In the diffusion term, VxSj = 



'^Sz^xSz comes from 

C(7i+V,7i+ + 7i_V,n_) = CV,(77^ + nl)/2 (13) 
and assumption of uniform n. The combination of trans- 



c- 



.dSz 
dx 



{2-qi + ■n2n)Ey = 0. 



(18) 



verse conductivity parameters rji and 7^2 in Eq.(ll) for 



It is easy to recognize that this equation describing the 
gradient of spin density reflects the balance of the diffu- 
sive spin current and the spin current generated by Ey 



This 



at 5*2 = 0, i.e., in accord with Eq.([5]), at a; = 0. 

is what defines the spin density, and, via Eq.(17), the 



transverse electric field in the whole central stripe, which 
is geometrically restricted by Eq.(14). Indeed, solving 
Eq.(18l with the boundary condition Eq. ^, we have 



2t]i+ ri2n 

Sz[X) = EyX 



(19) 



for the spin density. The spin density is the first order 
in Ey, and the first order in spin-orbit constants. The 
principal story is that the spin density grows or falls with 
X and inevitably reaches ±n. If it reaches n then n_ — 
O.If it reaches {—n) then n+ = 0. In both instances Eqs. 
( 14 1 are no longer valid. The coordinates —a < x < a at 



which this happens are the boundaries ±6 of the central 
stripe. 

For the transverse electric field in the central stripe we 
have 



E^x) - 



772(277: 



I^E'x 

2K '^ 



(20) 



We see that the electric field is quadratic in Ey and is 
independent of the polarity of the external field. Besides 
the skew scattering constant, it is defined by the com- 
bination of skew scattering and side jump contributions 
that describe the spin Hall current at a; = 0. The elec- 
tric field in the central stripe does not play an exclusive 
role, because transverse currents of electrons with each 
of the spin projections are also balanced by the diffu- 
sive spin currents of electrons with the corrresponding 
spin projection, due to spatially-dependent Sz{x) given 
by Eq.(f] 



The 



19) 



boundaries of the central stripe are given by 



± 



nC 



(2771 + rj2n)Ey ' 



(21) 



The boundaries ±& are defined by skew scatering and 
side jump currents via the combination that describes 
spin current at a; = 0. 

Once the maximal Sz = n ov minimal Sz = —n spin 
density is reached, it cannot grow (or decrease for op- 
posite x) any further, and Sz — ±?^ outside the central 
stripe, i.e., in the periphery stripes. Because the spin 
density is constant in each of the periphery stripes, and 
7T.+ = 0, n_ = 71 or n_ — 0, n_|_ = n, the diffusive terms 
in Eqs. ([2]) vanish. Spin currents in these two periphery 
stripes j± — a±{n±)Ey are compensated by the corre- 
sponding electric currents j = a{n±)Ex, and this gives 
two opposite electric fields in the two periphery stripes. 
The magnitude of this transverse electric field is given by 



\Ei'^\ 



(?7i -I- T]2n)Ey 



M 



(22) 



The combination of side-jump 771 term and skew scatter- 
ing 7^2 term here characterizes the magnitude of the AHE 
and spin Hall currents, which coinside at maximal spin 



polarization. On the formal grounds, the sign of the elec- 
tic field is opposite for the two periphery stripes due its 
dependence on the sign of the spin density, because Sz 



defines whether E[ 
3+ or i- : 



(p) 



is determined by the equation for 



sign{E^') — ~sign{Sz)- sign{rii+ri2ri)-sign{Ey). (23) 

We note that a convenient vantage point for the con- 
sideration of the electric field in the periphery stripes is 
afforded by Eq.(12). Indeed, because the diffusive term 



vanishes, it has a form 



lin 



Ei"^ + im + V2n)SzEy = 0, 



(24) 



and represents a balance of Ohm's current in the pres- 
ence of the transverse electric field and the AHE cur- 
rent that arises in periphery stripes due to the full spin 
polarization. At 15*2 1 = n, we arrive to Eq. (22 1 for 



the magnitude of the electric field. This picture looks in 
much the same way as in the conventional "charge" Hall 
effect^^, but here two opposite electric fields arise in the 
two periphery stripes. In order to find the direction of 

(v) 
Ex at positive and negative x, we determine the sign 

of the spin density on the boundaries of the central and 

periphery stripes, which determines the spin polarization 

in the periphery stripes: 

sign{Sz) — ''Sign{2r]i + 772^) • sign{Ey) ■ sign{x). (25) 

It is noteworthy that the sign of spin density in the pe- 
riphery stripes depends on the sign of the spin current 
at a; = (the centerline), governed by the combination 
(27;i +ri2n). Combining Eq.(23) and Eq.(25), we find the 

sign of {Ex ) defining the direction of the field in the 
periphery stripes: 



sign{Ex ) = [sign{Ey)]'^ ■ sign{'qi 
sign(2r]i + 77271) • sign(x). 



mn) 



(26) 



Remarkably, althouh the magnitude of the transverse 
electric field in the periphery stripes is linear in the ex- 
ternally applied electric field (and thus linear in the lon- 
gitudinal flowing electric current), the direction of Ex 
does not depend on the polarity of the external field. We 
underscore the cause of this effect: Ex' depends on Ey 
directly, and, independently of this, also on the sign of the 
spin density in the stripes. Having maximal (minimal) 
value, the spin density in the stripes does not depend 
on the magnitude of the external field but depends, as 
Eq.(25) shows, on the direction of the field. Thus, the 



sign(Ey) appears twice in Eq.(26), and upon the change 



Ap) 



of polarity of the external field, the direction of Ex does 
not change. However, the electric field is an odd function 
of the transverse coordinate x. 

The direction of the electric field in the periphery 
stripes turns out to depend on the product of the two 
combinations of side jump-like and skew scattering cur- 
rents. One combination defines the magnitude of AHE 



and spin Hall currents at maximal spin polarization, the 
other, coming from Sz on the boundaries between the pe- 
riphery and central stripes, defines the spin Hall current 
on the center-line of the sample. The transverse elec- 
tric field across the sample is defined by Eqs. (20 1,(22) 
and (26). A general feature of the transverse electric 



field is a steep dependence on coordinate at the bound- 
ary between the central and periphery stripes. In ze- 
roth approximation to charge-electric field distribution, 
at central-periphery stripe boundaries the electric field 
magnitude experiences "jumps" 



AS, = ± 



2771 



?72"- 



2/i 



E,, 



(27) 



The numerator in this equation is defined by the value 
of the spin current at x = 0. From the discussion of 
the electrostatic picture in the next subsection and Ap- 
pendix A, we will see that in fact there are two extremely 
narrow transitional ranges of x between the central and 
periphery stripes, in which the electric field has a finite 
derivative as a function of x. Still, if spin relaxation is 
negligible, from comparison of the value of electric field 
and its variation with x in the central stripe with the 
value of the electric field in the periphery stripes, one 
can separate the side jump-like currents from the skew 
scattering. 

In a real electronic transport experiment, a voltage 
rather than the electric field is measured. For uniform 
distribution of impurities, the electric potentials are the 
same at the two boundaries of the sample. However, 
voltage, which is directly related to spin current, can be 
measured between a probe at the edge of the sample and 
a probe at x in the periphery stripe. The magnitude of 
this voltage is 



f/.,n = 



(7?i + mn)Ey 



IJ- 



x\. 



(28) 



The magnitude of potential difference between the cen- 
terline and a probe at a point x within the central stripe 
is given by 



Uc{x) 



T]2{2rii +ri2n) 
4/xC 



'2^2 



E^x- 



(29) 



At x = the potential is chosen to be zero. This equation 
holds even if the sample is narrow, and the spin polariza- 
tion does not reach its maximal possible magnitude. The 
magnitude of potential difference between the centerline 
and a probe within a periphery stripe is given 



Upix) 



4/x(277i + T]2n) 



A* 



(30) 



The sign of these potential differences depends on electric 
field configurations. In particular, ii Ex {x < 0) < 0, the 
potential is higher at the edge than at the centerline. If 
Ex {x < 0) > the potential at the edge is lower than 
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FIG. 5: Dependence of spin polarization of electrons on co- 
ordinate transverse to passing electric current in case I. The 
slope of the dependence in the central stripe is calculated for 
parameters dicsussed in Sec. IV. Similar spatial dependences 
characterize cases II, III and V, but with the corresponding 
slopes and widths of the central stripes. 

that at the centerline. We discuss examples of these situ- 
ations below. If measured, potential differences between 
an edge of the sample and several probes within it would 
provide information about the magnitudes of side jump 
group currents or the skew scattering contribution. 

The overall shape of the transverse electric fields and 
potentials depend on the magnitudes and signs of the 
772 skew scattering and rji side-jump-like contributions. 
We consider four situations that will be described micro- 
scopically in Sec. IE with a justification of their physical 
parameters in Sec. IV. We mention here that the effect of 
iji currents has the same sign as the Drude conductivity, 
and the 772 current has the sign opposite to 771 contribu- 
tions in the case of domina nt scatte ring by the attractive 
impurity potential, see, e.g.,H31lMl. In the case of repul- 
sive impurity potential the sign of the 7/2 contribution is 
the same as the sign of the 7/1 currents. Then, in the case 
of remote doping with large setbacks, when the principal 
scattering mechanism is due to impurities in the quan- 
tum well, while the carriers are primarily due to remote 
donors, the skew scattering current can be eliminated by 
introducing full compensation of quantum well donors by 
acceptors. 

The qualitatively distinct cases are: 

• Case I. Skew scattering is absent: (772 — 0). The 
electric field Ex in the central stripe vanishes. In the 
periphery stripes, Ex^\x < 0) < and Ex^\x > 0) > 0. 
The spin polarization is p = 1 at a; < and p — —1 
at a; > 0. The spin polarization is shown in Fig. [5l the 
transverse electric field is shown in Fig. [6J the potential 
distribution is given by Fig. [2] 

• Case II. Side-jump-like term dominates both the 
AHE and the spin currents at a; = 0: 771 -|- 77277, > 0, 
27^1 -|- 7727T, > Skew scattering is characterized by 772 < 
0. The slope of the electric field in the central stripe 
dEx /dx > defines the camel back potential profile. In 
the periphery stripes, the situation is similar to case I: 
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FIG. 8: Electron potential energy profile along the direction 
transverse to the flowing current for case II. The underlying 
parameters are discussed in Sec. IV. 
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FIG. 6: Spatial dependence of transverse electric field for case 
I. The underlying parameters are discussed in Sec. IV. 
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FIG. 7: Dependence of transverse electric field on coordinate 
transverse to the fiowing current for case II. The underlying 
parameters are discussed in Sec. IV. 

Ei''\x < 0) < and Ei''\x > 0) > 0. The distribution 
of the spin polarization is also similar to case I: p = 1 at 
X > and p = — 1 at x < 0. In comparison with case I, 
the central stripe is wider and the periphery stripes are 
narrower. The transverse electric field is shown in Fig. [71 
and the potential profile is given by Fig. |8] 

• Case III. Side jump-like current is dominant both 
in the spin current at x = 0, 2r]i + 77271 > 0, and in the 
AHE, T]i+ri2'n > 0, and skew scattering contribution with 
772 > is present. The qualitative behavior of electric 
field and spin polarization in the periphery stripes is the 
same as in cases I, II. In the central stripe, the transverse 
electric field has a slope dEx /dx < 0, which results 
in a parabolic shape of the potential profile. The spin 
polarization across the sample is similar to that in Fig. [5] 
for the case I. The central stripe is narrower, and the 
periphery stripes are wider for the case III compared to 
cases I, II. The transverse electric field is shown in Fig. |9j 
and the potential profile is given by Fig. 10 



FIG. 9: Spatial dependence of transverse electric field for case 
III. The underlying parameters are discussed in Sec. IV. 

We will see in Sec. IV that cases I, II and III, describe 
the proposed experimental settings, although spin relax- 
ation modifies the spin density, electric field and potential 
energy shapes, as discussed in Sec. II. There are other 
shapes of spin density, electric field and potential that 
may characterize the spin-electric stripes in the absence 
of spin relaxation. It turns out, however, that in the 
InAlAs/InP/InAlAs system, which is a possible existing 
experimental setting for observation of the stripes, the 
ratio of skew scattering and side-jump-like contributions 
required for these shapes make spin relaxation strong. 
Still, we briefly discuss these shapes. 

• Case IV. Skew scattering 772 < dominates both 
the AHE and the spin currents at x = 0:771 -I- 77277. < 0, 
2771 + 77277 < 0. The profile of electric field and potential 
in this case is similar to case I. The spin polarization, 
however, reverses its sign compared to case I, as shown 
in Fig. [Tl] 

• Case V. Skew scattering (7^2 < 0) dominates the AHE 
771 -|- 77271 < 0, but the side jump-like contributions domi- 
nate the spin current at x = 0, 2771 -I- 77277 > 0. This is a 
surpprising situation. In the central stripe, the transverse 
electric field has a slope dEx /dx < 0. The spin polariza- 
tion is qualitatively the same as in the case I. However, 
the electric field in the periphery stripes is opposite com- 
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FIG. 10: Electron potential energy profile along the direction 
transverse to the flowing current for case III. The underlying 
parameters are discussed in Sec. IV. 
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FIG. 12: Qualitative picture of electric field profile along the 
direction transverse to the flowing current for case V. 
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FIG. 11: Qualitative picture of spin polarization profile along 
the direction transverse to the fiowing current for case IV. 



pared to the cases I-IV. The physical picture explaining 
this field distribution and the corresponding potential en- 
ergy shape is as follows. In the central stripe electrons 
with spin up travel towards negative x and reach the 
boundary of the periphery stripe with 5*2 = 1. In the pe- 
riphery stripe with these negative a;, the transverse cur- 
rent of particles with spin up, instead of continuing to 
flow towards the edge of the sample, as in Figs. (2[3 ), re- 
verses its direction and flows towards the boundary of the 
periphery and the central stripe. This is because the spin 
current at 5^ = 1 and the spin current at Sz = have op- 
posite signs for the set of parameters in the case V. Still, 
the sum of all contributions to the currents satisfies con- 
servation laws of Eqs. ( ll|12p in both the center-stripe 
and the periphery stripes. The corresponding profile of 
the electric field that compensates the reversed flow of 
spin (and AHE) current in the periphery stripes is shown 
in Fig. [12] and the potential profile is shown in Fig. [T3j 
Measurement of potential difference between the side 
contact to periphery region, and probes in the periphery 
and the central stripes is the most straightforward exper- 
iment. In this regard we note that the results presented 
above implicitly assume the dependence of the spin cur- 
rent on density given by Eq.( 10 ). While the term propor- 



FIG. 13: Qualitative picture of electron potential energy pro- 
file along the direction transverse to the fiowing current for 
case V. 

tional to 771 is universal (independent of impurity scatter- 
ing mechanism), we shall see that the skew scattering 772 
term generally depends on the ratio of the square of the 
transport time and the power 3/2 of the collision time 
(Sec. IE). For short-range scattering, when all scatter- 
ing times are energy-independent, the Eq.(lO) holds. In 



quantum wells with ionized impurity scattering mecha- 
nism being dominant, the skew scattering term in the ab- 
sence of screening could depend on density as n_(_ rather 
than rij. . Screening will bring back the power in the de- 
pendence closer to 2. A meaningful! experiment is then 
to measure the dependence of the transport time on den- 
sity, and use this dependence to solve Eqs. (15 16 ) with 
a modified density dependence of ri2 term. This will re- 
sult in modification of the solution in the central stripe, 
and a different functional dependence of the electric field 
and spin polarization there. The transverse electric field 
in the periphery stripes will remain linear in Ey, and will 



be defined precisely by Eqs. (22 1 with the only modifica- 
tion of density dependence of 772 term in brackets. It will 
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be then possible to separate the modified first and second 



terms in Eq.(30l and determine r]i and r]2- 



Other experiments aimed at observing spin-electric 
stripes are also feasible. For example, the electric field 
spatial distribution can potentially be measured optically 
via birefringence (more precisely, linear birefringence - 
the Pockels effeclP^. The spin structure also can allow 
testing via traditional measurement of spin polarization 
distribution. The periphery of the sample is 100% po- 
larized, the two domains have opposite spin polarization, 
and the polarization in each of the domains is constant 
in the absence of spin relaxation. We shall see in Section 
II that this picture is only slightly modified by the pres- 
ence of weak spin relaxation: the polarization reaches 
maximal possible magnitude only at the edge and de- 
creases within periphery stripes towards boundaries with 
the center stripe. We should also remark that only in ze- 
roth order in charge-electric field distribution, the total 
carrier density is constant: in the next approximation, 
discussed in the following subsection, the total density 
becomes non-uniform. Interestingly, while the circular 
polarization of the optical signal for the preiphery do- 
mains does not depend on the density distribution, the 
intensity of the signal can sense the modified charge den- 
sity. 

We underscore that maximal spin polarization caused 
by passing current through the sample is not limited to 
current-carrying electrons. All electrons in the periph- 
ery stripes are characterized by one z— projection of spin. 
The spin polarization at the centerline is zero, and is de- 
scribed by Eq.(fl9| in the central stripe. We reiterate 
that, strictly speaking, the above picture is valid when 
the z-component of charge carrier spins is conserved, so 
that spin relaxation plays no role in limiting accumula- 
tion of spins near the boundary, as well as in the balance 
of spin densities in the bulk. It is then the Coulomb 
interaction and the maximal possible spin polarization 
that set the limit on charge accumulation near bound- 
aries and leads to formation of the electric field stripes. 
Our solution for spatial distributions of electric fields and 
densities of each of the spin species takes into account the 
dependence of the electron energy and charge densities 
on the spin polarization. 



B. Electric field and charge density in the spin-electric 
state 

We now discuss the disctribution of electric fields and 
charge density using considerations of 2D electrostatics. 
So far, we have assumed carriers with both spins inde- 
pendent of each other and found that passing electric cur- 
rent results in transverse electric field periphery stripes 
with uniform and opposite electric fields. The electric 
field in the central stripe is linear in transverse coordi- 
nate, and depends on microscopic mechanism of the spin 
current. If only the side jump type currents are present 
then the electric field in the central stripe vanishes. We 



treat the electron density as constant across the sample, 
similarly to the conventional Hall effect. However, the 
electric field and the total charge density, which includes 
electron charge density n — n+ -I- n_ and a background 
charge density due to ions/impurities, define each other 
self-consistently. Indeed, the distribution of the electric 
field in the system is given by the Poisson equation 



diuE = 47rp, 



(31) 



where p is the charge density. Normal components of the 
field E2 are discontinuous at the 2D plane: 



S,(z + 0)-£;,(z-0) = 47r(iVo 



(32) 



where A^o is the 2D positive charge density (due to ions 
in metals or ionized impurities in semiconductors). The 
electron charge density n, depending on the experimental 
situation, comes either from electrons in metals, or from 
doping in semiconductors. In semiconductor hetcrostruc- 
tures, electrons come from doping of the 2D quantum well 
itself and from remote doping layers. 

We now demonstrate that in the next approximation 
to electric-field-charge distribution, the charge distribu- 
tion corresponding to the electric field stripe structure is 
non-uniform. For this demonstration, we consider the sit- 
uation when electrons come only from the quantum well 
impurities, so that the total charge of electrons and impu- 
rities in the 2D plane is zero (total electroneutrality) . To 
simplify our task, we will also restrict the self-consistent 
solution of the electrostatic problem to the case I with 
only side jump group spin currents present. We note 
that while the electric fields and the voltages between 
the center-region and the periphery of the sample are 
directly observable, the distribution of charge density is 
much more difficult to identify. True charge density de- 
pends on fluctuations in distribution of impurities. It 
would be of interest to address this and other aspects 
of the electrostatic problem, such as charge-field distri- 
bution in the presence of symmetric remote doping with 
two additional layers of positive charge, and to include all 
microscopic mechanisms of the spin current, for specific 
experimental configurations. 

We now solve Eq.( 31|32 | for case I, by using the bound- 
ary conditions for the electric field in the stripes and con- 
straints on electric charge density. The electric field in 
stripes is described by Eq.( 22|26 ) at 772 — 0, and is zero in 
the central stripe. The distribution of charge in a sample 
corresponds to the electric field obtained in the initial 
approximation to charge-electric field distribution, and 
is subject to constraint on the minimal value of possible 
electron charge. The origin of the constraint is as fol- 
lows. If the distribution of the electric field consistent 
with the spin current results in a discontinuity according 
to electrostatics, and the charge at discontinuity must 
be positive, then this charge density cannot exceed the 
density of positively charged ions (or ionized impurities 
located in a quantum well) A^o- We assume that each 
donor or ions supplies one electron, and correspondingly. 
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as soon as the Nq limit is reached, no electron charge 
can further be removed from a given location in order to 
reach values of positive charge bigger than Nq. This re- 
striction matters because, si milarly to the case of charge 
profile of the 2D charged stripSSEH^ -^ye can expect diver- 
gencies of charge density at the boundaries of the stripes, 
in which the electric field is non-zero. For example, in a 
2D charged metal film, divergencies of charge density at 
the edges result in a finite electric potential (f>. The pres- 
ence of excess electrons at the boundaries in our case is 
aided by the spin current. If electrostatics requires some 
ranges of x to be abandoned by electrons with both spin 
projections, there is a constraint on the remaining posi- 
tive charge at these x. 

In the long 2D strip configuration we only need to find 
the distribution of electric field in x— direction, and the 
dependence of charge density on x, because the system 
is uniform in the direction y of the electric current flow. 
Consider a complex function of complex variable E{t). 
The real part of this function Ex{x, z), ReE{x + iz) de- 
fines the electric field in the sample: 



2tt 



[E{x + iO) -f E{x - iO)] = E^{x), 



(33) 



and the complex part of —ImE{x -j- iz), as follows from 
Eq. ([32| , defines the distribution of charge: 



[E{x + iO)-E{x-iO)] = 4TT{No-n+{x)-n^{x)), (34) 

We now find the function E{t), using the following con- 
ditions: 

ReE{x) = Eo if h < |a:| < a (35) 



ReE{x) = 


if |a;| < b 


(36) 


Im[E{x + iO) - E{x - iO)] 


= -47rA^o 




ii b < \x\ < bi 




(37) 


Im[E{x + iO) - E{x - iO)] 


> -47rA^o 




if \x\ <b or bi <\x\ < a , 




(38) 



where 6i defines the range 6 < |a;| < &i of minimal pos- 
itive charge. These ranges are transitional regions be- 
tween the central and the periphery stripes. The value 
of 6i is defined from known Eq, Nq, by using the asymp- 
totics of E{t) as \t\ — > oo, which for electrically neutral 
sample gives E{t) —>■ 0(l/t^) as |i| — >■ c». The condition 
Eq. ( 37 1 is the constraint on the positive charge density. 
Furthermore, this constraint also affects the boundaries 



of ranges in Eqs (35J36J38). However, as we shall see, bi 
is exponentially close to b, and in reality the constraint 
is relevant only at the points b (positive and negative) 
merged with b-p^, which coinside with the boundaries 
between the central and the periphery stripes. By in- 
troducing Eqs.( 35|36|37|38 ) as above, we describe the 
electric field and charge distribution in the stripes self- 
consistently. 

The details of solution of Eqs.( 31|32 ) using methods of 
2D electrostatics are given in Appendix A. The electric 
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FIG. 14: Distribution of tlie transverse electric field due to 
the spin current across the 2D sample. A^'o = 5 x lO^^cjn"^, 
Eo = 2.4V/ cm. 




FIG. 15: Distribution of the charge accross the 2D sample, 
caused by the spin current. 2D charge carriers are assumed 
to originate from donors in the quantum well only, A^o = 
Eo = 2AV/cTn. 



field in the system is given by Eqs.( 163|164|167 ) and the 
charge density is defined by Eqs. ( 168|169|170 |. The 
electric field profile is shown in Fig. |14| (electroneutral 
quantum well). 

The charge density is shown in Fig. 



15 We note that 



the values b and &i that define intervals of x abandoned 
by carriers of both spin polarizations (this indeed occurs) 
are exponentially close, so that such intervals are essen- 
tially two singular points symmetric with respect to the 
centerline x = of the sample. However, in the vicinity 
of these points, there is a visible logarithmic dependence 
of charge density, which potentially can be detected. The 
logarithmic character of the position dependence of the 
charge density is defined by Eqs(T68jl70l, and Eq. ( 167) 
at rird = 0. We note that measurements of charge dis- 
tribution and spin polarization can potentially demon- 
strate the logarithmic behavior of density in the vicinity 
of |a;| = 6, and the square root features in the vicinity of 
the edges at |a;| — a. 



C. Spin-electric stripes in the presence of a perpendicular 
magnetic field 

Magnetic field H in z-direction (perpendicular to the 
2D plane) results in the two principal effects: (i) usual 
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(orbital) Hall effect resulting in a uniform electric field, 
and (ii) in a magnetic field induced spin polarization of 
charge carriers. This spin polarization, even if all other 
sources of spin polarization are ignored, results in the 
Anomalous Hall effect of spin-polarized carriers, i.e. the 
transverse charge current. In the constitutive equations 
([2]) for currents due to electrons with both spin projec- 
tions we must add the effect of charge Hall current, which 
at small magnetic fields is given by 



Jf 



IJ,n±9ujcT{z X E) 



(39) 



where Wc = eH/mc, is the cyclotron frequency. The coef- 
ficient 6 is defined by the energy dependence of momen- 
tum relaxation time r, and for r independent of energy 
6^1. In contrast to the spin Hall current, this contribu- 
tion comes with the same sign in equations for the cur- 
rents due to electrons with "-I-" and "— " 2;— projections 
of spins. Furthermore, in equations (pi), the densities n^ 
and n_ are now magnetic-field dependent. We note that 
taking into account this dependence of density on mag- 
netic field automatically accounts for the AHE current, 
and no further modifications of constituitive equations in 
magnetic field are necessary. 

In a degenerate 2D gas at low temperatures, in the 
absence of sources of polarization other than magnetic 
field, the electron spin polarization Sz — Pan, where 
Ph = gt^BH/ep-: g is the electron g- factor and /i^ is the 
Bohr magneton. Then on symmetry grounds, such spin 
polarization is the only source of the spin polarization at 
X = 0, and the boundary condition for Sz now is Sz(0) = 
Pnn. Solving the constituitive equations, for the spin 
polarization in the central stripe we get 



Sz\x) == pHn EyX. 

For the electric field in the central region we have 



(40) 



Ex{x) = -{pHfl- 



2ryi + ri2n 

c 



EyX)Ey+UJcTEy. (4l) 



The second term in this expression for the electric field 
comes from the charge Hall effect. It is interesting to 
note that similarly to the case H = 0, the side jump and 
skew scattering terms contribute to electric field in the 
central domain in a non-trivial way, which potentially 
gives a path for additional experiments in magnetic field 
that would allow their separation. 

We now consider the properties of the periphery stripes 
in the presence of magnetic field. The uniform Zeemann 
spin polarization results in the shifts of the boundaries of 
the stripes, and makes periphery domains asymmetrically 
positioned with respect to the centerline of the sample. 
Indeed, the boundaries ±6 between the periphery stripes 
and the central stripe are defined as x at which the elec- 
tron spin polarization reaches maximal possible values 
Sz/n ^ ±1. In the presence of magnetic field induced 
spin polarization, p^, depending on its sign, the bound- 
ary is closer to the centerline on one side and further from 
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FIG. 16: Dependence of transverse electric field on coordinate 
transverse to the flowing current in the presence of magnetic 
field perpendicular to the 2D plane. Case I. 

the centerline on the other side. The two boundaries are 
now given by 



(pH ± l)nC 
(2ryi + ri2n)Ey 



(42) 



Thus, application of magnetic field is capable of mov- 
ing the stripe boundaries. The electric fields in the two 
periphery stripes are given by 



Sf 



uJcTEy =F 



(771 -t- r]2n)Ey 



M 



(43) 



where the upper/lower sign corresponds to the spin polar- 
ization in stripes, Sz/n = ±1. As expected, the electric 
field due to charge Hall effect leads to electric fields in 
two stripes having different magnitudes. If the charge 
Hall effect field exceeds the magnitude of its spin Hall 
effect counterpart, electric fields in domains may become 
of the same signs, in contrast to their opposite signs at 
H — 0. For the charge Hall current having half value 
of the spin Hall current in the case I, the distriibution 
of electric fields in the sample is shown in Fig 16 The 



potential difference between the centerline a: = of the 
sample and the right edge x = a, U{0) — U{a) = Vh/^—V 
in the presence of magnetic field is different from the 
potential difference between x = and the left edge 
X = -a, U{0) - U{-a) = -Vh/2 - T^ by the amount 
of the charge Hall voltage, i.e. the potential differ- 
ence between the left and right edges U{—a) — U{a) = 
Vh ■ The half-sum of these two potential differences 
[2t/(0) - U[a) - U[-a)]/2 = -V defines the voffage drop 
between the edges and the centerline due to spin Hall 
current. 



D. Effects of non-uniform doping 

The ideal picture we considered so far is a metallic 
electrically neutral 2D plane. However, we shall see that 
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a feasible experimental setting is a semiconductor het- 
erostructure with remote doping. Remote doping leads 
to conditions of electrical neutrality fulfilled in the whole 
sample, so that the 2D plane itself may actually be 
charged. Were the 2D plane a metal, in the absence of an 
external current and therefore spin current, it would be 
an equipotential irrespective whether it is charged or not. 
Correspondingly, no fields in the Hall direction would ex- 
ist in the absence of external voltage, and, once the volt- 
age is applied, it would generate transverse electric field 
domains, in much the same way as we have discussed 
already. Thus, in the case of the charged metal the ob- 
servable voltage between the center of the sample and the 
edges clearly remains the same voltage shown in Fig. [5] 
Fig. [8] and [TO] The concentration profile in a charged 
metallic plane is going to be modified compared to the 
electrically neutral plane. We will not consider this issue 
here, and rather turn our attention to the semiconductor 
setting. 



In semiconductors, the quantum well is inevitably 
doped by impurities. We shall see below that in order 
to achieve the best conditions of experimental observa- 
tion of the stripes, it is useful to dope the quantum well 
itself, and in addition to use the remote doping. Dop- 
ing profile is never uniform, both in-plane and along the 
growth direction z. In particular, one of the effects of this 
non-uniformity (as well as non-uniformity of the remote 
doping), which we analyze below, is an additional chan- 
nel for the the spin relaxation of charge carriers. The 
question arises whether impurities can in any way alter 
the observable potentials (fields) in the proposed exper- 
iment, and obscure the observation of the electric field 
stripes. 



Non-uniform distribution of donors results in a built-in 
electric fields and the corresponding non-uniform electron 
charge density. Nevertheless, in the absence of an exter- 
nal voltage, and therefore absence of the spin current 
and transverse electric fields that it induces, measuring 
of voltage between the center of the sample and the edges 
results in zero signal in the steady state regime. The situ- 
ation is similar to measuring a voltage drop in a Schottky 
barrier or p-n junction. Without a battery in an external 
circuit, when the circuit is closed, one has two Schottky 
barriers (or p-n junctions) compensating each other, and 
zero signal. However, if the external voltage in longitudi- 
nal direction is switched on, the spin current is generated, 
and the transverse electric fields are induced in the pe- 
riphery stripes and the central stripe. Then the external 
circuit with a voltmeter will show the voltage drop pre- 
cizely as in Figs [2| [8| and [TOl Thus, while a non-uniform 
distribution of impurities in a semiconductor quantum 
well may alter the electrons charge density distribution 
compared to the ideal metallic case, the measured volt- 
age between a probe within the sample and a probe at 
the edge remains the same. 



E. Contributions to the spin current and their dependence 
on electron density. 

The microscopic picture of the spin current of electrons 
described by the Hamiltonian (fl]) includes two primary 
mechanisms: skew scattering assymetry and the effects 
associated with the change of the center of gravity of the 
wavepacket of electrons subject to electric field and im- 
purity potential (one of the effects in the latter group is 
the side jump current). It is important for us here that 
the sign of skew scattering current depends on the donor 
or acceptor character of impurities doping the quantum 
well, and that the magnitude of the skew scattering con- 
tribution is parametrically different from the magnitude 
of the effects in the side jump group, and has a different 
dependence on density of charge carriers. Skew scattering 
and side jump contributio ns ha ve been calculated in nu- 
merous works, see, e.g.,'^^!^^^'^. In order to explain con- 
ditions, in which spin-electric stripes arise or are absent, 
we need here the expressions for these currents in both 
3D and 2D conductors in certain experimental setups and 
spin relaxation times, in the same microscopic model. We 
therefore present the calculation of these quantities. 

The skew scattering contribution to the spin current 
density in the ideal case of Eq. (IT]) can be formulated in 
terms of distribution functions of electrons with spin up 
(or spin down), because there are no spin flip processses, 
and two spin species of electrons can be considered inde- 
pendently. We note that the effect of electron-electron in- 
teraction of skew scattering leading to the spin Coulomb 
drag can be also taken account by this approach. As 
was shown irF^, the spin Coulomb drug contribution ul- 
timately introduces a correction to the skew scattering 
contribution. Our estimates will show that the resulting 
reduction of the skew scattering effects amounts to just 
a few percent in the relevant range of parameters. 

The contribution to the spin current by skew scattering 
effect is given by 



Jsk 



e / d^kvk/r 



(44) 



where the spin-dependent (proportional to az) distribu- 
tion function J^" is obtained from the kinetic (rate) equa- 
tion for the electron distribution function /j, 



E 



dk 



k' 



iw^l 



+ W^k2')(/k-/k'), (45) 



where for 2D electrons 



<^ - 



27r 



|^k.k'P<5(ek-4), 



(46) 



Vk.k' is the Fourier component of the impurity potential, 
Ui is the impurity density, and S is the area of the 2D 
sample. For the calculation of longitudinal electric cur- 
rent and momentum relaxation time, it is sufficient to 
calculate the probability of scattering in the Born ap- 
proximation. 
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The probablity of asymmetric (skew) scattering in the 
framework of perturbation theory arises in the approxi- 
nation next to the Born approximation: 



W, 



(1) 



k.k' 



_ (2^Y 



X^k" -^"^[^kk"Vk"k'Vk'k] 



5{e^.~e'^)5{e^.~e'{,). 



(47) 



We underscore that the asymmetric skew scattering prob- 
abihty that we consider arises from generic 2D spin- 
orbit interactions given by Eq.Q and does not rely on 
the presence of intrinsic spin-orbit inter action s, e.g., of 
Rashba type, as in the case discussed iiPSEll. We will 
see that such contributions for various symmetries of in- 
trinsic spin-orbit interactions do not play any role in the 
context of spin-electric stripes induced by the nonequi- 
librium spin current in section IIIB.7. 

We will assume the model, in which the principal mech- 
anism of electron scattering is from point-like impurities 
doping the quantum well, 



y(r) - Y. ^^"^(^ 



R, 



(48) 



where Ri) gives the impurity location. The amplitude Vq 
is negative when Eq. ( 48 ) models attraction of electrons 
to ionized donors. The amplitude Vq is positive when 
Eq.(48l describes charged acceptors. Point-like scatter- 



ing is a reasonable approach if all charge carriers come 
from impurities in the quantum well itself. This implies 
that if charge carriers in 2D sample are supplied, apart 
from impurities doping the well, by the remote J-doping, 
then the remote doping layers should be separated from 
the 2D electron gas by sufficiently large setbacks. Then 
a smooth random potential from remote doping is negli- 
gibly small, and scattering by impurities of the quantum 
well are dominant. In this case from Eq. (fTl) the scatter- 
ing amplitude, including spin-dependent part, is 



Vkk' =V"o(l + iacr, -(kxkO^ 



(49) 



From Eq. (451, the asymmetric distribution function de- 



scribing skew scattering is given by 



k' 



>r{f, 



(E) 



I^') 



(50) 



where the momentum relaxation time r is calculated from 
Eqs.( 45|46|49 ) and in the leading approximation neglect- 
ing small spin-orbit term is defined by 



27rniVJ) V, 



(51) 



with 



2'Kh'^ 



being the electron density of states per 



spin; the distribution function f^ describing electron 
flow in a conductor in the presence of electric field is 
given by 



f{E) 
ik 



eAE.kr^Mfhl 
TO aek 



(52) 



where /o(ek) is the Fermi distribution function. Using 
Eqs.( 47|49|50|51 1, we obtain the expression for skew scat- 
tering contribution to the spin current 



Jsk 



-2an(j/3E, 



J" 



(53) 



where a = ne^T jra is the longitudinal conductivity, and 
/3 is the parameter distinguishing successive orders of the 
Born approximations: 



/3=^=± 



Txhri 



(54) 



where the upper sign corresponds to scattering by accep- 
tors and the lower by donors. The expression on the right 
of this equation can be used for estimates of /3, which can 
reach ^lHHH. 

When both donor and acceptor impurities are present, 
it is possible to suppress skew scattering at large compen- 
sation. Provided both donors and acceptors are shallow 
hydrogen-like impurities, at full compensation the skew 
scattering vanishes, but electrons coming from remote 
doping contribute to longitudinal mobility and side-jump 
group currents. Assuming the concentration of quantum 
well donors n^ and quantum well acceptors Ua, the total 
skew scattering current is given by 



J sk 



2ana(3KEy, 



(55) 



where K = (ud — na)/{nd + «-a) is the degree of com- 
pensation. Defined this way, K = Q corresponds to full 
compensation, and K = =pl correspond to only accep- 
tors or only donors present, correspondingly. We note 
that because of remote doping, the electron density in the 
quantum well n, and partial densities of electron species 
n± are the parameters of the system that can be tuned 
independently of K . 

Provided that r is independent of energy, the skew 
scattering current of carriers with a certain projection of 
spin polarization is quadratic in density of these carriers 
(being proportional to iJ±n±). The parameter 772 that 
determines the characteristics of the stripe state (Sec. 
lA) is given by 



m 



27r 



aPrK. 



(56) 



We now consider the parametric dependence of the ef- 
fects like side jump. These effects can be considered sep- 
arately for the two spin species of electrons in much the 
same way as the skew scattering contribution to the cur- 
rent. The side jump current is given by 



•>^.j=e^<,k'^k.k'/k, 



(57) 



k,k' 



where Wk,k' is the scattering probability (it is sufficient 
to take into account symmetric part Wq only), /k is the 
electron distribution function (it is necessary to account 
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for longitudinal flow of electrons in the presence of elec- 
tric field described by /^), and 



i?i 



k,k' 



.d_ 



■ 9 ^ 



Hi, Jii, 



(58) 



is the "side jump", i.e., the displacement of center of 
gravity of the electron wavepacket in the proceess of 
scattering^. Furthermore, $k,k' is the phase of the ma- 
trix element of scattering, and 



Hk = i 



d 



Ukg^UkdV 



(59) 



is the diagonal part of the matrix element of the elec- 
tron coordinate associated with the Bloch amplitudes Uk 
in the wavefunction. In the presence of spin-orbit in- 
teractions, Jlk = cio'z{kx,—ky,0) . This additional con- 
tribution to coordinate is often viewed as a dipole asso- 
ciated with the spin-orbit interaction, so that the spin- 
dependent part of the potential in the Hamiltonian (fl]) is 
considered to be due to this dipole moment. 

We would like to note here that the side-jump contri- 
bution expression (57 1 can be equaivalently derived by 
the three methods: (T) by taking into account terms that 
are second-order in scattering in off-diagonal density ma- 
trix and summing over all semiconductor banda^Zl (ii) by 
taking into account a single band and corrections of the 
first order in scattering to the current operato r toge ther 
with first order corrections to the density matrixSHSZI^ ^nd 
(iii) by considering a semiclassical picture of the contri- 
bution to the current arising from the shift of the center 
of gravity of the wavepacket in the course of electron 
scatteringSS. Apart from the direct contribution of the 



side jump to the current given by Eq.( 57 ), there are other 
similar contributions, which we collectively refer to as 
the side jump group effects. The side jump group effects 
remain under debate since 195423, and the common un- 
derstanding is that several contributions to those effects 
partially compensate one another, while the difference of 
opinions is to what extent this compensation occurs and 
what is the physical meaning of such compensation and of 
various contributions themselves. There is an agreement, 
however, what is the overall sum of side jump contribu- 
tions in the simple case of elastic scattering by impurities, 
as well as on the sign of the side jump contributions com- 
pared to skew scattering. Our goal here is to calculate 
the magnitude of the constant rji . 



Inserting (491, (51), (52) into Eq.(57l, we obtain 



Jsj 



ne^EyU 



- anE„ 



epT 



(60) 



As was shown irt^Z H5i | 69 ] ^j^jg ^^m-j^g Q^t ^q bg ^ correct ex- 
pression for the sum of all side-jump-group effects. Side 
jump current of species with certain spin projection is 
linear in density of these species as opposed to quadratic 
dependence on electron density of the skew scattering 
contribution. A signature of the difference between the 



two is the factor of ep in the denominator of the expres- 
sion describing jsj ■ The linear dependence on density 
for side jump currents holds even if the momentum re- 
laxation time depends on energy. The parameter rji that 
describes the characteristics of the stripes is given by: 



Vi 



2ea 
IT' 



(61) 



The side-jump-group currents can be made dominant 
over skew scattering contribution in a situation when 
compensation of impurities is close to 100% (i.e. almost 
equal number of donors and acceptors). While in 3D 
samples52J equal compensation is not possible because no 
current carrying electrons present in a sample, in 2D case 
the charge carriers can be supplied by the remote doping, 
and even full compensation of quantum well donors by 
quantum well acceptors is feasible. At certain compensa- 
tion, the skew scattering effect can vanish. If the strenth 
of the donor and acceptor potentials are the same, this 
happens at full compensation. As we show in Sec. Ill, in 
low mobility samples, the magnitude of spin relaxation 
rate is minimized, and that results in feasible experimen- 
tal setting. Results of calculation of contributions to the 
spin current are presented in Sec. IV. 



II. SPIN-ELECTRIC STRIPES IN THE PRESENCE OF 
SPIN RELAXATION 

Some spin relaxation will always be present in the sam- 
ple (even if Lso > 2a). In numerous contributions since 
the celebrated paper by M.I. Dyakonov and V.I. PereP, 
the effects of asymmetric skew scattering, and, particu- 
larly, spin Hall conductivity, have been discussed in terms 
of spin accumulation near the boundaries of the sample, 
which is limited by spin relaxation. 

The physical picture of electric field domains in the 
absence of spin relaxation is predicted in Sec. I. How- 
ever the question arises whether this effect can survive 
the presence of spin relaxation. I will now show that 
when spin relaxation is sufficiently strong, it is the spin 
relaxation that limits the accumulation of spins closer to 
the boundaries of the planar sample. However, at elec- 
tric fields larger than the certain field, the spin current, 
properly defined, reaches such a magnitude that the spin 
relaxation, and the associated spin diffusive current are 
not sufficient to balance the spin current induced y the 
external field, and limit spin and charge accumulation. 
Then the Coulomb interaction takes on the leading role 
and results in the transverse electric field and the spin- 
electric stripes. I will first demonstrate why spin relax- 
ation can be insufficient to limit the spin accumulation, 
so that electric field and the periphery stripes with 100% 
sample edge spin polarization must be a necessary in- 
gredient of the physical picture. I will then derive the 
properties of domains in the presence of spin relaxation. 
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Constituitive equations in the presence of spin 
relaxation. 



In the presence of spin relaxation, there is no longer 
a conservation of the currents caused by charge carriers 
with each spin projection, and instead, the spin density 
satisfies the continuity equation 



dx 



iUt - Jx ) 



2S, 



(62) 



A factor of two and the sign in the right hand side of 
Eq.(64| is due to our definitions of currents in Eq.(l2]) 
and spin density Sz in Eq.Q, and Tg is the microscopic 
spin relaxation time for the z— projection of spin. We 
assume the steady state and , for simplicity, no external 
magnetic field. The charge current is conserved: 



dx 



(j^ + j7) = o. 



(63) 



We will assume firs t that j^ and j~ are described by 
the same equations ( 2|6|9 10) as in section I.l, and thus 
consider the non- conservation of the spin currrent ( 62 ) 
to be the principal effect due to spin relaxation. This 
approach is along the lines ol^'. We will discuss the other 
possible roles of spin relaxation on a microscopic level in 
the end of this subsection. In order to pinpoint where 
the physical picture of spin accumulation gets replaced 
by the physical picture of spin-electric stripes, we will 
now use the same model for the source spin current as 
irP, i.e., assume linear in density spin Hall conductivity 

and set 7^2 = in Eq.( 10 1 . 

Inserting Eqs. ( 2|6|9|10 l into Eqs. ( 62|63 l, and assum- 
ing the total density n{x) — const and 772 — 0, we have 
the following equations for the transverse currents: 



_d_ 
dx 






2S^ 



0, 



(64) 
(65) 



The nonlinear terms in these equations were neglected 
in^, where no density or coordinate dependence of the 
coefficients in continuity equations were taken into ac- 
count. In this work, these factors are important. 



Wliy spin accumulation picture is insufficient 
and how spin-electric stripes arise 



The usual assumption for solving Eqs.( 64|65 ) is that 
Sz is proportional to the electric current, so that the 
first term in ( 64 ) becomes non- linear in electric field and 



therefore can be dropped. Let us follow this assumption 
for a moment. Then the same procedure applies to the 
S^ term, and all terms containing the electric field, and 
we arrive to the standard spin diffusion equation 



dx' 



;Sz 



Sz_ 

LV 



(66) 



where Li = (uts. This spin diffusion equation is solved 
with the boundary conditions at the edges x = ±a 



(^n d 
Ydx 



Sz + m^Ey = 0, 



(67) 



which reflects the absence of flow of the spin current out 
of the sample. We underscore that in the diffusion equa- 
tion itself, the gradient of the source term in this model 
vanishes, V xTjinEy = 0, but nevertheless the spin polar- 
ization arises due to the spin current affecting boundary 
conditions. 

What these equations describe is a traditional spin ac- 
cumulation picture. Let us now see whether this picture 
always applies. The solution of the spin diffusion equa- 
tion is 



(cosh{a/Ls) Lg 



(68) 



[To] 



This solution for a planar sample, see , e.g.Y^, is an ana- 
log of the 3D solution for a cylindrical sample irP. At 
small x (the center region) there is linear in x, Ey and 
771 spin polarization, due to the spin current. If the spin 
relaxation length is small, there is linear in electric field 
spin polarization at the edges, as very well known. 

However, let us now pose a question what happens if 
the spin relaxation length and time are sufficiently long, 
and the magnitude of electric field is sufficiently large, to 
the extent that 



\Sz{b,\b\<a)\^±n, 



where 



L^arsh 



'Cncosh (a/Lg)^ 
■. rjiEyLg 



(69) 



(70) 



This amounts to the spin polarization reaching its max- 
imal possible value at some x such that < a; < a (and 
minimal possible value for the opposite spin projection at 
> —X > —a). Then at 6 < |a;| < a, the solution given 



by Eq.(68) not only results in spin polarization exceeding 



its maximal possible value, but in contradiction to this 
value being maximal possible, the solution continues to 
grow with X approaching the boundaries. For negative 
polarization, at the opposite coordinate, Sz/n < —1, and 
the spin polarization according to the solution decreases 
with X approaching the boundary despite having already 
reached the minimal possible value. Therefore, the so- 



lution ( 68 1 fails at sufficiently long spin relaxation times 



and large electric fields. 

The reason is, in these conditions one cannot neglect 
the electric field in Eqs.( 65|64 ). Spin relaxation of the 
z-component of spin alone cannot limit the accumulation 
of spin caused by the spin current, and the electric field 
takes over this role. We are going to show now that the 
range of coordinates —b<x<b defines the central stripe 
and the ranges —a < x < —b and b < x < a define the 
spin-electric periphery stripes. 
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C. The solution of constituitive equations for 
spin-electric domains 

To find the elecric field and the spin density, we must 
solve the system of the non- linear equations ( 64|65 ) sup- 
plemented by the boundary conditions of the absence of 
flow of the spin current and charge current out of the 
sample at a: = ±a: 



fJ-^z^x H Y'dx^ 



finEx 



4 dx'^z 



+ VinEy 

riiSzEy -- 







(71) 
(72) 



We will assume that spin relaxation is weak. In the model 
we discuss now, the corresponding condition is 



T]iEyn > 



Cn2 
L. 



The more general condition at 772 7^ is 
-1-^ D+n+ 



(73) 



(74) 



The physics of this condition is that the spin current 
caused by the electric field is much stronger than the 
diffusive spin current caused by non-uniformity of spin 
density on the scale of spin relaxation length. We im- 
mediately note that the possibility to reach this regime 
is non-trivial. Indeed, if to assume that the spin current 
and spin relaxation are governed by the same spin-orbit 
constant a, then 771 oc a, 1/ L^ ex a, and the strength 
of spin interaction does not enter this condition. If that 
were true, too strong electric fields, in which electrons 
gain energy eEyt ^ ep on the scale of the mean free path 
£, are required to fulfill (73). However, we shall see that 



in 2D gas different spin constants govern the spin current 
and spin relaxation. 

We now will solve the system of equations ( 64|65 ) by 
perturbation theory using the small parameter of Eq. 
(73 1. In the zeroth approximation, i.e. l/r^ — 0, the 
stripe solution of subsection I. A. holds. We next find 
first order corrections in I/t^ to this solution. In the pe- 
riphery region, say, around x = —a we write down the 
spin polarization in the periphery stripe in the form 



S^ix) ^n- s{x), 
and the electric field in the form 



Ex{x) = - 



VlEy 



+ e{x) = Eo + e{x), 



(75) 



(76) 



where s{x) and e{x) are first order in l/r^, and Eq is 
the electric field in the stripe at the boundary x — ~a 
in the absence of spin relaxation. Inserting Eqs.( 75[76 ) 
into Eqs.( 64|65 ), and keeping only linear in 1/t^ terms, 
we have 



V 



, ( - fiEos{x) + Tnie{x) - ^ V,s(.x)) = "^ 
VxUme{x) - ^Vxs{x) - riis{x)Ey\ = 0. (77) 



This system must be solved with boundary conditions on 
spin and charge current for e{x) and s{x), which amounts 
to the two expressions in round brackets in ( 77 ) vanishing 



at a; = —a. This turns out to be equaivalent to 
s{-a) = 0. 



(78) 



We note that because both the electric field and spin po- 
larization are the odd functions of x, it is sufficient to 
solve the system of equations at x < in order to ob- 
tain the whole solution. At negative x for the periphery 
region, we have 



s{x) 



n{x-\-a) 

~ VlEyTs ' 



e{x) = 
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The resulting spin polarization and the electric field in 
the periphery stripe are given by 



SP ^n I 



E. 
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(81) 
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where b is the x-coordinate of the boundary of the central 
and the periphery stripe in the absence of spin relaxation, 
given by Eq.(21) at 772 = 0. The absolute value |5| is the 



half-width of the central stripe, and at cc < 6 is nega- 
tive. Compared to the picture with spin relaxation ab- 
sent, the magnitude of the electric field is slightly smaller 
at a; = —a, and it decreases further as x arroaches the 
center region, as opposed to being constant. The spin 
polarization has its maximal possible value at the edge, 
and instead of keeping this value throughout the periph- 
ery stripe, decreases towards the center region. 

We now find the solution of Eqs. ( 64|65 ) in the center 
region, in the linear in l/r^ approximation. For 7/2 = 0, 
we have found E^ = in the absence of spin relaxation, 
as seen from Eqs.( 17|20 |. We assume that the electric 
field can arise in the first order in l/r^: 



Ex{x) = ec(x). 
For the spin polarization, we assume 

2mEy 



Sl{x) 



c 



Sc{x), 



(83) 



(84) 



where the first term is the spin polarization in the cen- 
ter region in the absence of spin relaxation, and Sc{x) is 
the first order correction in spin relaxation rate. Then, 
inserting Eqs.( 83|84 ) into Eqs. ( 64|65 ), we obtain the 
following equations for edx) and Sc{x): 

V. ( - ^-^x ■ e,{x) + 'f V,s,(x)) = -^a;(85) 

fj-nedx) - riiEyVxScix) = 0. (86) 

These equations must be solved with boundary condi- 
tions reflecting that the spin polarization and the electric 
field vanish at x = 0: 



Sc(0) 
gc(0) 



0, 

:0. 



(87) 
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The solution of these equations for Sc{x) is : 



Sc{x) 



1 nx 



(89) 



and therefore the spin polarization in the center region is 



^^^(-) = -^-(l + S)- 



LI 



The solution for the electric field is 



E^{x) 



2Eobx 



(90) 



(91) 



Here the product of Eq and b is positive for in the whole 
range of x. We note that this result (as our whole pro- 
cedure) applies in the range of electric fields compati- 
ble with Eq.(73) , so that the limit i?j, — >■ automat- 



ically implies I/ts — )■ 0, resulting in zero electric field 
in the central stripe. Formally the result (91) originates 



from the non-linear term Sz^xSz in the equation for the 
charge current, which includes the product of the two 
contributions in spin polarization given by Eq.(90|. The 



non-linear term itself follows from the dependence of dif- 
fusion coefficients for electrons with two spin projections 



on their corresponding densities, as defined by Eq.(13) 



Next, from the continuity of spin polarization we find 
the boundary bs of the central and the periphery stripes: 



5^(5,) = S^,{b.). 
Using Eqs.( 81|90 l we obtain 



.-.n-m^.-%-' 



(92) 
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In the limit of strong electric fields b^ and b coinside, 
because b oc 1/Ey. We note that the second bracket 
formally gives divergence of bg when the width of the 
central stripe is equal to spin relaxation length. This is 
precizely the condition when our expansion at small spin 
relaxation rate becomes invalid, as expected. 

For the set of parameters discussed in Sec VI, the spin 
polarization, the electric field and the potential energy 
in the presence of spin relaxation are shown in Fig. |19[ 
Fig. [20] and Fig. [21] correspondingly. We observe that the 
stripe state manifests itself indeed even in the presence 
of spin relaxation. There are three distinct regions, the 
center stripe and the two periphery stripes, with spin 
polarization reaching maximal possible value (positive 
or negative) at the edges of the sample. In the initial 
approximation to charge-field distribution, the value of 
the electric field in the presence of spin relaxation ex- 
periences jumps from the central stripe to the periphery 
stripes. Applying considerations of electrostatics as in 
Sec. IB and Appendix A, we find that the derivative of 
the electric field is not infinite (albeit large) as it changes 
from the center to periphery. The corresponding field 



small correstions linear in I/ts- Similarly to the case of 
absence of spin relaxation, the change occurs in a very 
narrow range of x close to bg . The boundary of the cen- 
tral stripe and the periphery stripes at finite Ts is marked 
by an abrupt change of the VxSz{x), as it occurs in the 
absence of spin relaxation. 

In our consideration of the stripe structure in the pres- 
ence of spin relaxation, we introduced several simplifica- 
tions. Like in the case of absence of spin relaxation, we 
assumed that there are two independent channels of con- 
ductivity, diffusive currents and spin currents, and the 
only effect of spin relaxation taken into account has been 
the modification of "continuity" equation for spin cur- 
rent via 1/ts term. In general, there are contributions to 
conductivity and diffusion, which are off-diagonal in spin, 
caused by spin relaxation, but the leading terms in dif- 
fusion and conductivity are those we took into account. 
These off-diagonal contributions can change the details 
of corrections to spin polarization and electric field in 
the stripes, but not the principal features. Furthermore, 
spin relaxation can result in additional contributions to 
spin currents (or charge currents of spin-polarized carri- 
ers). For str ong sp in relaxation these contributions can 
be significantpSl,CSI. We analyse several possible contri- 
bution to spin and AHE currents in system under con- 
sideration in Sec. IIIB.7. The upshot is that tn the 
case we consider, the condition for existence of stripes 



Eq.(73) requires that the constant responsible for spin 



current is significantly bigger than the constant leading 
to spin relaxation. It then follows that contributions to 
spin current due to spin relaxation in the relevant regime 
will always be small. Therefore our approach is justified. 
In our illustration of domains in the presence of spin re- 
laxation we discussed how the picture of domains caused 
by spin current defined by rji changes when I/ts is fi- 
nite. This allowed direct comparison with traditional 
spin accumulation picture. If we take Sec. lA stripes 
in the absence of spin relaxation defined by both rji and 
772 (side-jump-like and skew scattering terms), then, al- 
though the details of change itroduced by finite I/t^ will 
differ, the principal conclusion is still that in spin relax- 
ation gisve only a minor perturbance of the stripe picture. 
The next important task that needs to be accomplished 
is to prove that it is indeed possible to make the spin 
relaxation small but the spin current big. 



III. STRENGTH OF SPIN RELAXATION IN 

VARIOUS ELECTRONIC SYSTEMS AND 

FEASIBILITY OF ELECTRIC FIELD DOMAINS 

A. Spin relaxation and spin current in the 3D bullc 

systems 

1. Spin relaxation in the 3D bulk systems 



In bulk systems, the condition Eq.(74) requires very 



in the leading approximation is given by Eq.(162| with strong applied electric fields that are comparable with 
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or larger than internal confining fields of a heterostruc- 
ture in semiconductors, or, in metals, give energies ac- 
quired by electrons in electric field on the scale of mean 
free path significantly exceeding the Fermi energy, which 
would lead to melting of the sample. As we are going 
to show, the reason is that the 3D spin-orbit interaction 
Hamiltonian in a metallic film 



•H'" ^aa-[kx Vrl^(r)] 



(94) 



includes the spin-flipping terms containing Ca;, Cj,, which 
result in the relaxation of the z-projection of the electron 
spin, and this relaxation is defined by the same spin-orbit 
constant a that defines the spin current. 

We first discuss the effect of spin relaxation caused 
by spin fiip in 3D samples. General expressions applied 
for 3D case will also be needed for quasi-two-dimensional 
samples. The latter are relevant for for the spin relax- 
ation in quantum wells, where it is important to take into 
account a random distribution of impurities along the 
growth direction of the potential well. If spin-dependent 
scattering becomes a factor in spin relaxation in the case 
when usually more dominant effects in 2D samples are 
suppressed, it is this effect of the thickness of the quan- 
tum well that leads to the existence of the spin-fiip scat- 
tering, in contrast to the ideal 2D case described by Eq. 
([l]). To address these questions rigorously, we consider 
the evolution of the 2x2 spin density matrix p in the pres- 
ence of spin-dependent scattering. In 3D samples, the 
spin relaxation is defined by the Hamiltonian ( 94 1 . The 



quantum kinetic (rate) equation describing the evolution 
of the spin density matrix pk in the presence of spin- 
flip scattering can be derived by the standard Keldysh 
diagrammatic technique, see, e.g.,'25Ell^ a^d is given by 
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dp 
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(95) 



where [A, B] = {AB + BA)/2. Scattering of electrons 
with average spin projection Si (i.e., electrons described 
by the component of the density matrix proportional to 
Siai with Si independent of the direction o of k) in the 
presence of spin-orbit interactions ( 94 1 results in the two 



effects: (i) spin relaxation of electrons described by spin 
relaxation time Tso, and (ii) the appearence of the correla- 
tion between the electron spin and momentum described 
by the spherical harmonic (1 — 3of). The effect (ii) is 
described by the component of the density matrix pro- 
portional to Siai{l — 3o|). Here we are interested in the 
spin relaxation only. The spin relaxation time for i-th 
projection of spin given by 



1 



'^so,k'. 



All 



E (i^kVP " i^kV.p)^(i^k - i?k')), 



^ (96) 
where "H^^/ ^ is the part of matrix element ( 94 1 propor- 
tional to the i-th component of the Pauli operator, and 



angular brackets denote averaging over the directions of 
o. The spin relaxation time Tso is isotropic in directions 
of o, and, in general, is energy-dependent. For point-like 
impurity scattering, with scattering amplitude defined by 
Eq. ( 48 ) and modified Eq. ( 49 ) that includes terms propor- 
tional to (Ti • [k X k']j, i ^ x,y, z, we have 
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(97) 



where the scattering time r for 3D case is determined by 
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r3D 



V^n, 
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Then the spin-orbit Icgnth Lgo 

J^so ^ V J-JTso ^^ 
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(98) 



(99) 



where vpT = I is the mean free path. As clear from 
Eqs.( 97[99 ), the spin relaxation rate for this mechanism 
grows as a square of the Fermi energy and is proportional 
to the rate of electron-impurity collisions. In bulk semi- 
conductors, this mechanism is called Elliott- Yafet mech- 
anism, while in metals the analogius effect was discussed 
by Overhauser as early aa^S. In the case of 3D metals 
with inversion symmetry, Eqs.(98 99) fully define spin 
relaxation properties. If the symmetry of the crystall 
is lower, in heterostructures or in confined geometries, 
the Dyakonov-Perel mechanism of spin relaxation in 
the presence of intrinsic spin-orbit interactions is often 
the dominant spin relaxation mechanism in semiconduc- 
tors. In terms of electron spin resonance line width, 
the spin relaxation rate given by Eqs.( 97|98 ), character- 
izes the effects of collisional broadening of the line. The 
Dyakonov-Perel spin relaxation mechanism is an exam- 
ple of motional narrowing effecta^S. What is important, 
the collisional broadening spin relaxation mechanism in a 
bulk 3D metal cannot be suppressed once electrons con- 
tribute to the conductivity and scatter of impurities. In 
Sec.IIIA.2 we will see that this renders spin-electric do- 
mains (which are bulk analogs of spin-electric stripes) 
impossible. However, in quasi 2D samples, the situation 
is different, and we shall see that spin relaxation caused 
by spin-orbit scattering can be suppressed. We shall also 
see that intrinsic spin-orbit coupling due to low symme- 
try and the Dyakonov-Perel spin relaxation can be nearly 
eliminated. In order to discuss the spin-flip scattering in 
the 2D case, we need to derive the corresponding Hamil- 
tonian. The corresponding spin relaxation time will be 
discussed in Sec. IIIB.5. 



2. Spin current and external electric fields required for 
appearance of domains in bulk 3D systems 

We now calculate the spin current for bulk systems. 



Using ( 94 ) , for the contribution of skew scattering effect 
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we obtain 
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(100) 



where a refers to contribution of one projectrion of spin 
to conductivity. The side jump contribution is given by 
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(101) 



We note that for 3D metals, when the carrier density 
considerably exceeds the concentration of impurities, the 
square root in ( 100 ) is of order of unity, and at tpT /fi 3> 1 
currents of side jump group are much smaller than skew 
scattering contribution. 

Using Eqs (74 100 99), we calculate the electric field 



that is needed to generate electric field domains induced 
by spin current in the bulk system. In order to discuss 
feasibility of an experimental observation of the effect, 
we multiply this electric field by the mean free path and 
the electron charge: 



eEvpT 



mVf 



(102) 



With the denominator in this expression reaching mini- 
mum when the square root « 1,^3 it follows that energy 
acquired by charge carriers in electric field on the scale 
of the mean free path exceeds the Fermi energy. These 
are unrealistic conditions, certainly impossible in met- 
als. Note that if the side jump group contributions in 
some metallic conductors (metallic conductivity implic- 
itly means epr/h > 1) exceed skew scattering conductiv- 
ity, the electric fields needed for domains result in ener- 
gies acquired by electrons on the scale of the mean free 
path is epr/h times bigger than the Fermi energy. The 
current induced by such external fields would melt the 
sample. Thus, the spin relaxation due to ordinary bulk 
spin orbit interactions effectively prohibits appearence of 
the electric domains caused by the spin current. 

We observe that in order for electric fields to be experi- 
mentally feasible, energies acquired by charge carriers on 
the scale of the mean free path needs to be parametrically 
smaller than the Fermi energy. We now demonstrate that 
it is possible to find quasi-2D systems, in which applied 
electric field needed for observation of this effect is rea- 
sonable. These systems are characterized by a sizable 
spin current but the spin relaxation there is eliminated 
altogether, or at least suppressed by several orders of 
magnitude. 



[110] crystallographic direction. Spin relaxation in these 
systems has been already experimentally shown to be 
much smaller than for structures grown along [OOlp^ESl. 
Let us now discuss if there is the ultimate limit for spin 
relaxation rate in such setting. 



1. Intrinsic spin-orbit effects that need to be suppressed 

Spin-orbit interactions in quantum well structures, and 
semiconductors with low symmetry, apart from spin-orbit 
scattering given by Eqs.( l|94 |, generally includes intrin- 
sic spin-orbit effects, i.e. terms in the Hamiltonian 
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-?ia.(fl«+n(3)). 



(103) 
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Here n^' = Y^^Nkp^^, ^ = T.,utl^3Mk,kuku 
where i,j,u,t denote projections on the principal axis of 
the crystal x \\ [100], y |i [010] and z \\ [001]. In bulk 
III-V systems, cubic in k terms are present with ^xxyy — 

"^XXZZ — ^yyzz — ^yyxx — ^zzxx — ^zzyy and all 

other 7 components vanishing. This is the Dresselhaus 
terrrPS. Linear in k terms can arise in bulk crystals such 
with wurtzite or tellurium symmetry' 'l Linear in k 
terms in bulk III-V systems arise only in th e pre sence of 
uniaxial strain described by strain tensor ei P ' [ 



mntr = ^.-<T = 




Ci{axkx{eyy ~ ezz) + (Jyky{ezz - 


- e.xx) 


'^zkz{exx - f^yy)) + C2(exy{o'xky - 


— (JyKx 



+<^yzieTykz - (T,^ky) + Czxic^zkx " CT^^/c^)) (104) 

We note that non-zero diagonal components of the 
strain tensor naturally arise in the context of lattice- 
mismatched heterostructures. In the absence of strain, 
linear in k terms in III-V systems, such as GaAs, arise 
in 2D samples. In particular, the bulk Dresselhaus terms 
upon quantization along [001] direction to 2D become 
2D Dresselhaus terms with /3^ = — /3™ (all other f3^ 
are zero). In quantum wells with x,y conducting plane 
perpendicular to growth direction z, subject to asymmet- 
ric potential confinement there arises Rashba term with 
l^xy = — /3j,J^. AH these terms potentially result in spin 
relaxation of the z-component of spin. As we shall see 
below, in 110 grown III-V crystals, the situation is differ- 
ent, and intrinsic spin-orbit interactions leading to such 
spin relaxation are almost non-existent. 



Intrinsic spin-orbit interactions for 110 orientation of 
the quantum well 



B. Constructing quantum wells with no z-spin 
projection relaxation due to spin-orbit interactions. 

The example of the setting, in which it is feasible to 
observe electric field domains is the n-type semiconduc- 
tor double heterostructures (quantum wells) grown along 



As we have seen above, rather moderate spin relax- 
ation precludes spin-electric domains or stripes due to 
spin current. Our task is therefore to find systems with 
spin-orbit interactions such that the terms responsible 
for spin current are sufficiently strong, but the terms re- 
sponsible for spin relaxation of the z-projection of spin 
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are suppressed or vanish altogether. Let us consider het- 
erostructures grown in [110] direction. In this case, our 
z'-axis is along [110] crystallographic direction, a;'-axis is 
along [110] crystallographic direction, and y' \\ [001]. In 
these systems, several spin-orbit effects vanish**^. In par- 
ticular, the bulk Dresselhaus cubic in momentum terirp2I 
results only in the terms proportional to the combination 
of spin operators describing the spin projection along the 
normal to the 2D gas z': 
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(kx'K , k^i/z). 



(105) 



The Dresselhaus 2D spin-orbit interaction arising from 
averaging over the quantized motion alon g the z'- 
direction of the bulk cubic Dresselhaus ternPlMl also 
gives only terms containing spin operator projection 
along z': for electrons of the ground level of size quanti- 
zation, assuming for simplicity an infinitely deep rectan- 
gular quantum well of width d, 
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These terms give neither spin relaxation of the z'- 
projection of spin, nor the spin currenlPZl. If the quan- 
tum wells have symmetric heterostructure confinement, 
another source of spin-orbit interactions, the Rashba 
temPSl vanishes. Furthermore, if neigboring heterostruc- 
ture layers have common atoms at the interface, yet an- 
other source of the intrinsic spin-o rbit in teractions, the 
interface asymmetry, does not aris#2I^. Moreover, for 
lattice-matched systems, a strain does not result in ex- 
tra spin-flipping interactiona^Sl as well. Indeed, similar to 
the case of Dresselhaus terms in [110] grown structures, 
strain-induced terms Eq.(104| in the presence of the rel- 
evant strain tensor component ez'z' contain only a single 
Pauli matrix <7z' '■ 
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These terms also do not result in spin relaxation, and 
do not induce the spin current. Thus, the only possible 
source of spin-orbit effects in this system is the impurity- 
induced or fluctuation-induced spin-orbit interactions. 



3. Impurity and fluctuation-induced spin-orbit interactions 

These interactions are of two types. Interactions of the 
first type, described by Eq.([I]), conserve the spin projec- 
tion perpendicular to the 2D layer. However, because of 
the finite thickness of the 2D layer, there is a possibility of 
a second type of interactions, where the impurity-induced 
terms lead to spin relaxation. Impurities are randomly 
located with respect to centerplane z = of the quantum 
well, and the potential gradient of the impurity potential 
at the impurity location is sizable. Then, in much the 
same way as in the 3D conductors, at the first glance the 



spin-flip scattering is feasible, and can lead to spin relax- 
ation. Furthermore, it has been suggested by Sherman^ 
that random distribution of impurities in doping layers, 
even if it is symmetric with respect to the centerline of 
the quantum well, results in some minimal spin-orbit in- 
teraction of a "random Rashba" type that leads to spin 
relaxation. 

In the Appendices B and C we derive the effective spin- 
orbit Hamiltonian in systems confined to two dimensions 
in the presence of 3D random potential. The dependence 
on the third dimension turns out to be important. Nev- 
ertheless, we demonstrate that the effect of spin-orbit 
scattering by impurities or the spin-orbit effects in the 
presence of the long-range random potential on electrons 
in the quantum well can be greatly suppressed compared 
to similar effects in the 3D systems. Moreover, we show 
that these sources of spin relaxation can be eliminated 
altogether. 

The results of appendices B and C are as follows: In 
the presence of 3D potential V{x,y,z) affecting conduc- 
tion electrons*"*, there are three principal contributions 
to the spin-orbit interactions. 

• Spin-conserving spin-orbit scattering, i.e. the ideal 2D 
spin-orbit Eq.(fl]). 

• Spin flipping spin-orbit scattering governed by the val- 
ues of band offsets 

• Spin flipping spin-orbit scattering governed by the dif- 
ference of semiconductor parameters in the quantum well 
and the barriers. 

The division of parameters that control spin-orbint in- 
teractions and the related relaxation of the z-projection 
of electron spin in quantum wells into two groups is a 
matter of convenience. The first group includes band 
offsets in the valence and the conduction band; the sec- 
ond group includes all the parameters that could dif- 
fer in the barriers and the quantum well: the effective 
mass, the Kane parameter of the Kane model for III- 
V semiconductorgSSl (and potentially parameters of band 
coupling in other multiband models), and the spin-orbit 
splitting of the valence band. 

a. Band offsets and disorder-induced spin-orbit inter- 
actions In this paragraph we formulate the problem and 
present the results for the band-offset controlled spin- 
orbit constants. Details of this consideration, includ- 
ing construction of the effective 2D Hamiltonian starting 
from 3D interactions, and the role of boundary conditions 
are discussed in Appendix B. 

We derive a specific expression for spin-orbit constants 
for electrons in symmetric quantum wells that are of pri- 
mary interest in this work. Our general conclusions apply 
to arbitrary confinement. For simplicity we use the no- 
tation x, y, z for Cartesian coordinates, with z being the 
growth axis. We keep in mind the symmetry restrictions 
coming from the [110] growth direction of the quantum 
well structure. However, our results can be applied to 
spin phenomena in quantum structures with an arbitrary 
crystallographic orientation. 

Electrons in the conduction band of the double het- 
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FIG. 17: Potential profile of the heterostructure of the second 
type, in which electrons are quantized to a quantum well, but 
holes are not. Different features of the potential in barriers 
(subscript b) and the quantum well (w) are shown. 

erostructure are affected by confinement and scattering 
potentials 

U{z) = lJc{l~e{z+d/2)+e{z^d/2)) + V{x,y,z), (108) 

where d is the width of the well, Uc is the conduction 
band offset. The valence electrons are subject to the 
potential (whose relevance is discussed in Appendix B) 

C/„(z) = -E'^ + U^{\-e{z+d/2)+e{z-d/2)) + V{x,y,z). 

(109) 
Here E"]!^ is the band gap in the quantum well, and Uy 
is the valence band offset. We note that the only re- 
quirement is that conduction electrons are confined to 
form a 2D setting, while valence electrons for positive 
offset can experience a barrier potential and no confine- 
ment. This latter setting realizes a so called type If 
heterostructure^^ . 

For direct gap semiconductor III-V structures, the 
structure of valence band includes two bands: Fg band 
of heav y and light holes, whose confinement is described 
by Eq. ( 109 1, and the split-off Fy valence band, in which 



the potential acting on valence electrons is described by 
Uv7 = Uv — ^(z), where the coordinate-dependent spin- 
orbit splitting is given by 
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where A;,, A„ are energy separations of the split-off band 
in the barrier and the quantum well, correspondingly. 
The heterostructure potential profile is shown in Fig. [TT] 

Taking into account the difference of offsets in the con- 
dution and Fg valence band, and therefore the difference 
of the bandgaps in the barrier and the quantum well 
(Eg' — Eg ^ Uy ~ Uc), but neglecting the difference of 
all other parameters in the well and the barriers, we ob- 
tain the following effective Hamiltonian (including both 
spin-conserving and spin-fiipping terms) 
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where {A)oo = [/_oo dz'i'* {z)A{z)'i! {z) is the expectation 
value of the operator A{r) in the electron ground state in 
the confining potential U{z). The constant a describes 
generic spin-orbit interactions in the Kane model, and is 
given by Eq. ( 181 1 . With our choice of ordering of oper- 
ators in Eq.(112), a is negative. 



The first term in Eq.(112) has the same symmetry as 
spin-orbit interaction (HI) in an ideal 2D system. The 
two other terms provide channels for relaxation of the 
z— component of spin. Assuming wavefunctions ^(2) to 
be real, the spin flipping terms take a simple form 
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aK[{F{x, y)[pxal + [px a],F(x, y)]/2, (113) 



where the parameter k = 1 — Uy/Uc, and the scatter- 
ing potential gradient averaged over the ground state of 
the quantum weU is F{x,y) = (VzT^'"^°*)gQ. The sym- 
metrization here gives a Hermitian Hamiltonian, arising 
from an accurate averaging over the growth direction 
(Appendix B). The term coming from the differential op- 
erator kj^ = ~id/dr± acting on F{x,y) is especially im- 
portant in the case of scattering potential U{x, y, z) due 
to quantum well impurities. 

We now make several observations about these spin- 
orbit interactions. The spin-flip scattering terms in 



Eq.(112) vanish if potential offsets in the conduction 



band and in the valence band are equal. We underscore 
that in majority of heterostructures offsets differ signif- 
icantly and have opposite signs. Then k enhances the 
spin-flip scattering rate, as opposed to suppressing it at 
close offsets. 

We now discuss the origin of the dependence of the co- 
efficient K in spin-fiip spin-orbit scattering on potential 
offsets. Part of this spin-orbit term comes from aver- 
aging of dUy/dz over the ground state in the quantum 
well. If this averaging would only include averaging of 
the dVscat/dz over the wavefunctions calculated in the 
absence of scattering, then k = 1. However, scatter- 
ing potential disturbes the electron wavefunctions in the 
quantum well. Then it is necessary to include the con- 
tribution of the gradient of symmetric confinement po- 
tential averaged over scattering-disturbed wavefunctions. 
That changes k. Because the gradient of symmetric po- 
tential in rectangular quantum well is defined by offset 
potentials, k depends on offsets. Furthermore, a spe- 
cific combination of offsets in k comes from the structure 
of perturbation expansion in the Kane model, in which 
electron spin-orbit effects are related to admixture of the 
valence band potential. Moreover, because Vgcat is essen- 
tially a part of the overall conflning potential, vanishing 
K for Uy = Uc is related to the Ehrenfest theorem. In- 
deed, according this theorem, the average gradient of the 
conflning potential vanishes in any conflned state. For 
equal offsets, the spin-orbit term under consideration is 
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exactly such an average. We note that because Hamilto- 
nian describing spin-orbit interaction must be Hermitian, 
its part containing kzdUy/dr± is also proportional to k. 
In discussion of the 2D spin-related effects, such as 
weak antilocalization, or 2D spin relaxation in scientific 
literature, the term containing —idF{x,y)/dr± has not 
played any role. This has been a consequence of a com- 
mon neglect of the dependence of scattering potential on 
the coordinate normal to the 2D gas, giving a vanish- 
ing diagonal matrix element of Pz in the ground state 
in Eq.(112). Also, it has been usually asumed that the 



barriers and the quantum well. These contributions to 
spin-orbit interaction that arise due to 3D asymmetric 



Dyakonov-Perel spin relaxation due to the Dresselhaus 
or Rashba term is dominant. However, when intrinsic 
spin-orbit interactions are significantly suppressed like in 
the growth axis [110] configuration, the scattering terms 
are important. 

b. Spin- orbit interactions due to a difference of pa- 
rameters in the quantum well and barriers We now dis- 
cuss the contributions to spin-orbit interactions associ- 
ated with differences in materials parameters in the bar- 
riers and the quantum well. As follows from Appendix C, 
the spin-orbit interaction Eq.(112| arising in the model 



of constant effective electron mass in the double het- 
erostructure has to be modified if we take into account 
that simultaneously with potential offsets, the effect of 
variation of the value of materials parameters takes place 
across the structure. In particular, there are different 
masses in the barrier, ms, and in the quantum well, raw- 
For different masses in the barriers and the well, one ob- 



vious modification of ( 112 1 is simply to take into account 



the effect of mass difference on the wavefunction of the 
confined state ^(z). However, the principal result is that 
the additional contribution to spin-orbit interaction due 
to mass difference arises, which is entirely independent 
from Eq.(112). This contribution does not contain the 



parameter k that vanishes at equal offsets, but instead 
contains the factor ruB — rnw : 

(iV,vI/*^(d/2)[k X a]z^Vz^.,y{d/2) - 
iV.vI/*^M/2)[kx fT].iV,vI/,^M/2)), (114) 

where fluxes —Vz'^xyiz) are conserved across single het- 



erojunctions. The expression in brackets in Eq.(114) 
is nonzero because the scattering potential V{x, y, z) or 
asymmetric potential UasynAz) admixes odd and even 



states of the Hamiltonian (197), and makes ^2,y(z) a lin- 
ear combination of these odd and even states. We note 
that that the wavefunction '^xy{z) has the dimensionality 
l/-\/d, so that SV-^so has the dimensionality of energy. 

In Appendices B and C we have also considered the 
variation of the Kane model parameter P and spin-orbit 
splitting A between the Fg band of heavy and light holes 
and the split-off Fy across the double heterojunction, and 
derived the corresponding contributions to spin-orbit in- 
teractions caused by difference of these parameters in the 



potential, and scattering effects are given by Eq.(216|. 

c. Remarks on Rashba interaction due to asymmetry 
of quantum well confinement We now discuss the con- 
sequences of our consideration for Rashba effect due to 
asymmetric potential confinement or applied bias. Al- 
though in order to observe spin-electric stripes one has 
to get rid of asymmetry of confining potential and de- 
sign the doping layers symmetrically with respect to the 
quantum well, so that intrinsic spin-orbit interactions of 
the Rashba origin are removed, our results are useful for 
consideration of other phenomena due to spin-orbit in- 
teractions. When only Uasym{z) is present, in the ab- 
sence of random impurity potential, the Eq.(114| sim- 



plifies, so that different electron masses in the barriers 
and the quantum well result in an additional term of the 
Rashba symmetry 



5Ul^ = ah^niB - mw)"^^ 
(iV,**(L/2)iV,v|/(L/2)- 

iV,vI/*(-L/2)i-V,vI/(-L/2)) [k X a]. 



(115) 



where the ground state wavefunction of an electron in 
asymmetric potential '^{z) is independent of x and y 
coordinates. Similarly, for other contributions to spin- 
flipping spin orbit terms discussed in Appendix C, such 
as terms due to different spin-orbit splittings of the va- 
lence band or different Kane parameters in the barri- 
ers and the well, for ID asymmetric potential along the 
growth axis the Eq.(216) simplifies and becomes a term 
of the Rash ba symmetry given by Eq.(221). Indeed, all 
of our Eqs. pT2l fTTal fTwl [2161) describe not only the ef- 



fect of random potential V{x, y, z) on electrons confined 
to a symmetric heterostructure by the potential U{z) 
but also the case when there is weak asymmetry of con 
fining potential or externally applied voltage, leading to 
Rashba terms. The Rashba term arising fr om E q 
was discussed ir#^. The additional terms (1151, 
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221 



|) 

found here in the presence of asymmetric potential lead 
to a significant contribution to the Rashba interaction. 

We note that in a recent work by Chalaev and 
Vignale^, Rashba spin-orbit interactions arising from 
Coulomb interactions have been discussed. The authors 
assumed that confinement potentials acting on conduc- 
tion electrons and valence electrons are proportional to 
each other, and deduced that Coulomb interactions is 
the only source of the Rashba term. While in general 
this conclusion is in line with our result that scattering of 
electrons by impurities or by potential fluctuations can be 
the dominant source of spin-flipping terms, we would like 
to offer a few clarifications to^^. First, it was assumed 
that confinement potentials for valence and conduction 
electrons are proportional to each other (AlGaAs/GaAs 
heterostructure was considered as an example). This is 
certainly true for the offsets part of potentials in double 
heterostructure quantum wells. However, for the asym- 



metric potential, whether it is related to the external gate 
voltage, or to some smooth variation of the bandgap, this 
is generally not the case. Therefore, sources of asymme- 
try of the overall potential acting on valence electrons, 
be it an arbitrary position of an impurity with respect 
to the centerplane z = of the quantum well, exter- 
nal voltage, asymmetric confinement, asymmetry in the 
location of delta-doping layers, or fluctuations of long- 
range potential for symmetric doping, will all result in 
the spin-flip terms. These effects will still all vanish once 
the difference of materials parameters in the barrier and 
the quantum well is ignored. Incidentally we note that 
consideratiorPSl did not include the dependence of the 
Coulomb potential on the in-plane coordinates, which 
might be of interest in certain effects due to spin-orbit 
interactions, such as spin dephasing in quantum dots^^. 

d. Feasibility of spin-orbit interactions with a forbid- 
den flip of the z-projection of spin It is important to 
note, that while simple Rashba symmetry terms can be 
removed by making the confinement potential symmet- 
ric, the terms caused by the spin-orbit scattering due to 
impurities and monolayer fluctuations derived here are 
always present when the bandgap or other materials pa- 
rameters change across the heterostructure. 

These contributions can be significant especially when 
materials in the barrier and the well differ substantially 
in their Eg, A and P. Thus, it is clear that for suppress- 
ing the spin relaxation one needs to choose a structure, 
in which the parameters in the well and the barrier(s), 
particularly Eg, are somewhat close. This latter condi- 
tion is in line with the condition that n almost vanishes. 
Indeed, if one engineers a heterostructure in which Eg is 
constant, the spin relaxation would occur because of a 
minor possible difference of the Kane parameters in the 
barriers and the quantum well, while the spin-flipping 
terms in Eq.(fl2| would vanish precisely. Such confine- 



ment of electrons (in the absence of hole confinement) 
can be achieved in particular in systems with no abrupt 
hetero-boundaries but rather in situations when the con- 
fining potential for electrons, and a gradually varying po- 
tential in the valence band that repells holes are identical. 
In this case, a smooth variation of the band edges iden- 
tical in the valence and the conduction bands, as shown 
in Fig 18 leads only to existence of spin-orbit terms (II]) 
but not to spin-flipping terma^Sl. This is a particular ex- 
ample of realizing symmetry that gives spin current but 
no z-projection spin flips. 



4- Spin-orbit interactions due to large-scale potential 
fluctuations. 

Two types of disorder effects can contribute to spin- 
orbit interactions ( |112[ ). A random potential V{x,y,z) 
comes from not only a short range disorder and charged 
impurities in the quantum well, but also from a long- 
range disorder, particularly large-scale fluctuations of im- 
purity density resulting in a smooth random potential. 




FIG. 18: Ideal quantum structure confining conduction elec- 
trons but not holes with constant bandap and identical pa- 
rameters in the well and barriers, in which the efltective Hamil- 
tonian of spin-orbit interactions is given by Eq. (II]) 

We now consider spin-orbit effects associated with the 
large-scale fluctuations. 

The large scale fluctuations were discussed in^^ for 3D 
spin-orbit effects and in^^l in the context of minimal spin- 
orbit interactions in 2D gas confined by a symmetric po- 
tential. Our Eq.(112|, adds the following important fea- 
tures to this spin-orbit effect (i) the factor k = 1 — U^/Uc 
appears as a consequence of the rigorous procedure for 
averaging over the ground state in the present work; (ii) 
the third term of ( |112 [) which makes the Hamiltonian 
Hermitian is missed in^-. An accurate averaging over 
the z— direction is important because we need here a 
precize cap on spin-orbit interaction strength and spin- 
relaxation. We show that such averaging gives additional 
suppression of spin-orbit effects due to long-range poten- 
tial fluctuations. 

We consider the spin-orbit interactions in the quantum 
well located in the range of coordinates —d/2 < z < 
d/2 caused by impurity Coulomb potentials of the two 
narrow [5) doping layers of width w located in the range 
of coordinates —I — w/2 < z < —l-\- w/2 and —I — w/2 < 
z < —I -\- w/2. The 3D ionised impurity density in the 
layers is 7ii(R), i — 1,2. Average over the realizations 
of impurity potential, denoted by angular brackets, gives 
(rti(R)) — (n2(R)) = (n), so that the 2D density n2D m 
the quantum well induced by the doping layers is n2D — 
2w{n{'R.)). We assume the correlation function of the 
impurity densities in the doping layers to be of the "white 
noise" type: 

(in,{R) - (n))(n,;(R') - (n)) = (n>J(R- R') 

((n2(R)-(n))(ni(R')-(n)) = (116) 

The 2— component of the electric field due to impurity 
density fluctuations acting on electron on the first level of 
size quantization (with coordinate r = {x,y,z)) is given 
by 



6E,{x,y)^J2j'^^J^^^'^^( 

i 

( 



^^dr^^I^)^''^^^^"^''^^- ^^^^^ 



This fluctuation-induced electric field results in a ran- 



dom Rashba term defined by Eq.(113) with F{x,y) 
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6Ez{x,y). The average field Ez{x,y) itself vanishes be- 
cause of the syminetrically located doping layers. 

The correlation function of the fluctuating electric 
fields {SE{x,y)SE{Q,0)) is expressed in terms of aver- 
age square fluctuation electric field {5E{0, 0)^) and the 
correlator C(x,y)221 : 



{SE,ix,y)SE,iO,0)) = {SE,{0,Of)C{x,y). 



(118) 



Assuming here for simplicity the quantum well with infi- 
nite potential confinement with the wavefunction '^(z) — 



j cos ^, for {SE,{0,Oy) we obtair 



{SE,{0,Or)^ 



Trn2De'^l3^ 



(1 - 2l/3[ci{2lf3) sin (21(3) - si{2lf3) cos (21/3), (119) 

where /3 — 27r/d. At large AiTl/d, i.e., for large setbacks, 
this expression reduces to 



{SEMOr) = 



7rn2De 
2Pe2 



(120) 



Thus, a precise averaging procedure over growth direc- 
tion reduces {SE{Q,0)'^) compared to^'^. For the correla- 
tion function C{x,y), at setbacks large compared to the 
width of the quantum well, we obtain 



C{x,y) = 



1 



iir 



3/2- 



(121) 



where p = \/ x'^ + y^ . 

If the average Rashba interactions are suppressed by 
symmetric configuration of the two doping layers that 
are sufficiently close to the quantum well, non-uniform 
fluctuation-induced Rashba interactions may become a 
dominant spin-orbit coupling in high mobility materials, 
with the possibility to realize experimental settings sug- 
gested irPS. Here we are interested in suppressing these 
interactions, which is possible via large setbacks of dop- 
ing layers. 



5. Spin-orbit interactions in 2D systems. Suppression and 

enhancement of spin relaxation due to spm-orbit 

interactions. 

Once several intrinsic spin-orbit interactions in the 
structure are eliminated, the spin-orbit scattering in a 
quantum well with the scattering rate given by Eq.(96) 



becomes an important z- projection of spin relaxation 
mechanism. For spin-orbit interactions defined by 
Eq. ( |112[ ) , the matrix element of spin-flip scattering due 
to short-range interactions with impurities (48 1 in the 
ground state in the quantum well is given by 

((fc, + fc;)<'-' - (fc. + K)<j-y-'y^i<^-^'). (122) 



Here we included for brevity only the spin-orbit interac- 
tion related to offsets, deflned by k = 1 — ^. Spin-orbit 
terms considered in Appendix C associated with a dif- 
ference of parameters in the quantum well and barriers 
result in a contribution to the total spin-orbit constant 
that is added io an. 

We note that it is accurate inclusion of both spin- 
fipping terms in Eqs.( 112|113 ) that makes this scattering 
amplitude Hermitian, 



(«"):: - ii^'O 



(123) 



Inserting the matrix element (122| into Eq.(96), and av- 



eraging over impurity locations across the quantum well 
Z — i, we obtain the spin relaxation time 



1 



2d2 



(124) 



The corresponding spin-orbit length is given by Lgo — 
{dkp){ppT / m) / apjpK. 

The presence of the factor k, or factors describing small 
variation of semiconductor parameters from barriers to 
the quantum well provides an important method for 
technological control of the strength of scattering-related 
spin-orbit effects, average Rashba spin-orbit interactions, 
and fluctuation-induced spin-orbit interactions. 

Let us now discuss whether structures with strongly 
suppressed spin relaxation readily exist. Indeed, there 
is a system which is very close to the ideal case. In 
Ino,52AlQ,4aAs/InP/Ino,52Alo.4sAs double heterostruc- 
ture, for which we obtain n ~ 0.33 using known band 
offs ets an d band-gaps in the barriers and the quantum 
.yyg} ]94 | 95 | Such value of k decreases z— projection spin 
relaxation rate associated with (112) order of magni- 
tude (assuming kp '^ 1/d) compared to the spin relax- 
ation of transverse (x,y) spin components. Furhermore, 
the effective electron mass in InP differs from that in 
I no, 52 AIq AS As by just 5%, which makes the contribution 
of Eq.(114) to spin-orbit interactions very small. Other 



heterostruture systems can be hopefully engineered with 
even smaller k. We note, however, that in several other 
III-V compounds, z— projection spin relaxation in 2D 
case can be stronger than the relaxation of the transverse 
components of spin, due to the parameter k, or because 
of a sizable difference between the quantum well and the 
barrier parameters. In particular, the z— projection spin 



relaxation rate associated with terms ( 122 1 is the order of 
magnitude stronger than that for x— and y— components 
in 6*00.47/710.55^5 quantum well. 

We note that in GaAs/AlGaAs system, where a con- 
siderable decrease in the spin relaxation rate has been 
observed recentljBSl in the case of the Dresselhaus and 
average Rashba terms of equal magnitudes, the random 
Rashba contribution and spin-orbit scattering induced re- 
laxation of the z— projection of spin cannot be completely 
eliminated using the mechanism of reduction of k dis- 
cussed here. This is because the experimental setting^Sl 
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implies a sizable Rashba term equal to the Dresselhaus 
term., and thus nonvanishing k. However, as follows from 
Eq.(|120|), random Rashba term due to fluctuations can 



always be suppressed by increasing a setback (the dis- 
tance between the quantum well and each of the two 
doping layers equally separated from it.) 



6. The Dyakonov-Perel spin relaxation for non-uniform 
intrinsic spin-orbit interactions 

The spin-dependent terms that contain Ux and Uy re- 
sult in spin relaxation of the z-component of the spin. If 
such terms are independent of the electron coordinate or 
their coordinate dependence is smooth, the probability 
of a spin filp in a single scattering event is small, and the 
effective mechanism of the spin relaxation takes its origin 
in the Dyakonov-Perel spin diffusion^. Terms in ( 103 1 , 



( 104 1 constitute an effective momentum-dependent field 



that makes electron spin precess. In the absence of scat- 
tering, this precession gives two spin-orbit "chiral" quan- 
tum states defined by on the projection of the electron 
spin on the effective field. Ordinary spin-independent 
scattering due to random impurity potential Uir) leads 
to a change of the electron wave vector, and thus ran- 
domly changes the direction of the effective momentum- 
dependent magnetic field. If the angle of electron spin 
precession between two consecutive scattering events is 
small, the resultant effect is a random spin walk, dif- 
fusion, leading to spin relaxation. For the coordinate- 
independent effective field, this phenomenon is very well 
discussed in the literature from different vantage points. 
In particular, it can be described in terms of a phase 
breaking of the electron wavefunction due to effective 
spin-dependent vector potential, which makes transpar- 
ent the role of intrinsic spin-orbit interactions in weak 
localization and related phenomengPSHSll. 

In the context of spin-electric stripes, when sev- 
eral spin-orbit effects are suppressed, the coordinate- 
dependent intrinsic spin-orbit interactions can play an 
important role. In order to compare the strength of var- 
ious contributions to spin relaxation, we present an ac- 
count of Dyakonov-Perel spin relaxation taking this co- 
ordinate dependence into account. 

In the presence of intrisic spin-orbit interactions, and 
spin-independent impurity scattering, the electron spin 
relaxation is described by the quantum kinetic equation 
for the spin density matrix 






'^intr(l")'/'k} 
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^i5(ek -ek')|C^kk'P(/3k -Pk' 



Rashba interaction. We describe this interaction in terms 
of the effective spin-orbit vector potential Aso{v±) : 



ftintr 



pAso(r_L) + Aso(r_L)p 
2m 



(126) 



where for brevity we dropped a subscript c in our 
notation for the effective mass m, y±^ is the in- 
plane coordinate of the quantum well, A.so{yi_) = 
a2mSEz{x,y){ay,ax,0), oli = an -\- 5a, where 5a is the 
spin-orbit constant defined by differences of parameters 
of materials in the quantum well and barriers, Ez{x,y) 
is described by the Eqs.( |117|118|120|121[ ). We note that 
the coordinate- and spin-dependent vector potential also 
describes the spin relaxation in mesoscopic ferromagnets 
with domain walls and other magnetic textures^^. 

In the relaxation time approximation for a collision in- 
tegral due to impurity momentum scattering, the quan- 
tum kinetic equation ( 125 1 describing electron spin re- 



laxation and spin diffusion can be written as 



|-fP.VrP-P[A..(r),p(r)]=-^ 
at m ra r 



Po 



(127) 



where square brackets give the commutator of operators 
that depend on Pauli matrices, and po is the momentum- 
independent component of the density matrix. Here we 
drop the index _L for coordinates. 

As was shown in^°°, the electron-electron scatter- 
ing, which does not contribute to relaxation of mo- 
menta in semiconductors (that happens only in metals 
due to umklapp processes^, contributes to the effec- 
tive relaxation time that defines spin relaxation of elec- 
trons via the Dyakonov-Perel mechanism. The effect 
of electron-electron scattering also affects the electron 
spin relaxation for non-uniform intrinsic spin-orbit in- 
teractions. Therefore r in Eq.(127) includes a contri- 



bution from the electron-electron scattering defined by 
T'^ — T^ +T~J-. For degenerate semiconductors at low 
temperatures t^,} = C2ei?/?i(r/ei?)^, with C2 = 3.4^221, 
Defining Fourier-components p(r) = J dqp{q) exp iq • r, 
Aso(i") — J c?qAso(q) cxp iq • r, we obtain 



dqi[Aso(q- 



gp(q) 

dt 

qi),p(qi)] = 



— p • qp(q) - 

m 

p(q) - p(q)o 



(128) 



(125) 



In order to describe spin relaxation, we consider an ex- 
pansion of the spin-dependent part of p(q) in harmonics 
of p. In the presence of Rashba-like interactions, it is 
sufficient to restrict this expansion to terms independent 
of Pi and terms linear in pf. 



In order to consider spin relaxation caused by residual 
intrinsic spin-orbit interactions due to long-range fiuctu- 
ations of impurity doping, we now discuss the Dyakonov- 
Perel spin relaxation due to the position-dependent 



p{q)^a-S^"\q) + a,Sl]\q)p, 



(129) 



Then, inserting Eq.(129) into Eq.(128), multiplying by 
az/2 and taking the trace over spins, we express Sz (q) 
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in terms of S, 



(1). 



9S(°)(c 
dt 



(1), 



qSl^'p, 



-tr^ /dqi^[A,„(q-qi),a.5^;^(qi)p,] (130) 



(1)^ 



:.(i) 



(0) 



In turn, 5^ (q) is expressed in terms of S^ (q) 



asL"(q) 



9iC 



(1)^ 



at ' m'-'ij ^'i' 

-TmJrfqii[A,„,,(q-qi),a,5i°'(qi)] 



(131) 



In the semiclassical limit, the term containing a com- 
mutator of the effective spin-dependent vector potential 
Ago and the spin components describes the spin pre- 
session with frequency ft — a2kpEz{v)/h. We assume 
that the classical angle of rotation is small, Or <C 1 dur- 
ing the time interval equal to the momentum relaxation 
time. Using this condition and inserting the solution of 
Eq.(131) into Eq.(130), we obtain the equation of the 



i(0)/ 



evolution the average spin Sz (q)- In order to obtain 
the magnitude of the spin relaxation time, it is sufficient 
to examine the resulting equation in the stationary state, 
which reads 



2c0 



Dq^S' 



S° 



where 



= alkj,{6Ez{0,0f 



(132) 



(133) 



D = VpT/d is the diffusion coefficient, and d is the di- 
mensionality of the system. In the model that we use 
for a calculation of the spin relaxation caused by the 
fluctuations-induced Rashba term, the spin diffusion co- 
efficient D turns out to be the same as the electron diffu- 
sion coefficient. This is because t is defined by the spin- 
independent scattering, and therefore, S^^'pj relaxes in 
the same way as the electron momentum. Once the spin- 
dependence of scattering is taken into account, this will 
no longer be the case. However, in systems with potential 
spin-independent scattering being the dominant scatter- 
ing mechanism, the difference between the spin diffusion 
coefficient and the diffusion coefficient is negligible. 



It is worth noticing that as seen from Eq.(133) for 



the fluctuation-induced Rashba term, the spin relaxation 
rate is expressed in terms of average squared of the elec- 
tric field caused by fluctuations, (^£^^(0,0)^), and does 
not depend on the correlation function C{x, y). 

In much the same way we obtain the spin relaxation 
time in cases when impurity scattering depends on scat- 
tering angle. The procedure is also easily takes into 
account the smoothly varying intinsic spin-orbit inter- 
actions of various symmetries. If \^i\, \^i\ given by 
Eqs.( 103|104 | are much less than the frequency of impu- 
rity scattering, the solution of equation for spin density 



matrix in the second order in the effective precession field 
gives the spin relaxation time 



''s 



(134) 



((f2(3)2 + f^(3)2))r3, 

where -^ = X]k'(^^0)k,k'(l ~ cosnS), 9 is the angle 
between k and k'. For point-like scattering, ti = rs = r. 



7. Analysis of additional contributions to the spin and 
anomalous Hall currents 



The spin and AHE currents discussed in Sec. IE are 
due to spin-conserving spin-orbit interaction given by 
Eq.([T]). In our discussion of systems with reduced re- 
laxation of the z— component of spin, we have seen that 
several other spin-orbit terms are present. These are: (i) 
spin-flip scattering described by Eq.(122|, (ii) the ran- 
dom Rashba interaction given by Eqs. ( 113|126 ), and 
(iii) intrinsic spin-orbit interactions described by Eqs. 
( 106|105|107 ) containing only z— Pauli operator. The 
question arises whether these effects can result in a siz- 
able contribution to the transverse current. 

We begin with the discussion of the spin-flip scatter- 
ing. Using Eq.(122|, it is easy to see that the spin-flip 



processes result in both skew scattering [Eq.(53)] and 
side-jump- like [(Eq.(57l] currents. However, these con- 



tributions are quadratic in the spin-orbit constant re- 
sponsible for spin relaxation, an. An example of the 
structure of such contribution, e.g., for skew scattering 
of electrons with a spin projection a, is as follows: one 
of the matrix elements in Eq.(47), which now includes 



the summation over spins of intermediate states, is pro- 
portional to kx{(JyY''^ , where the spin projection a' is 
opposi te t o a. Of the remaining two matrix elements 
in Eq.(47l one includes a term k'{axY ''^i ^nd the other 
is spin-independent. The resultant sum over a' contains 
kxk'y{(JzY''^ ■ Such a term leads to generation of the spin 
current in x— direction in the presence of non-equilibrium 
distribution function /^ ex k' Ey. Because this effect is 
quadratic in q;k, its ratio to the transverse current due 
to spin-conserving interaction Eq.([I]) is described by the 
parameter 

^ ~ 0.7 X 10"^ (135) 

where the numerical estimate is made for the In- 
AlAs/InP/InAlAs quantum well with 300A width. As 
we seek a system with weak relaxation of the z— com- 
ponent of spin, every extra order of spin orbit constant 
results in additional smallness of the corresponding con- 
tribution. Thus, we can neglect these effects. 

Discussing the role of random Rashba interactions, we 
first note that for uniform linear in momentum spin- 
orbit interactions, no intrinsic contributions to non- 
equilibrium spin conductivity is generated. This fact was 
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derive d in n umerous contributions over the past decade, 
see e.g. OS— "231, as well as is apparent from the calculation 
of the component of the nonequilibrium spin density ma- 
trix proportional to the Pauli operator and momentum 
in early work^^ . Because of smooth non-uniformity of the 
random Rashba term, a small spin current arises, which 
at weak spin-orbit interactions is quadratic in an. There- 



fore, its smallness is characterized by Eq.( 135 1, and it can 
be neglected. Second, random Rashba interactions can 
serve as a sourse of asymmetry that leads to skew scat- 
tering or to side-jump-like contributions, instead of spin- 
orbit term in Eq.(ll]). Such contribution was discussed, 
e.g., for skew scattering AHE current iiiSasfiij fgj. gtrong 
uniform Rashba interaction, and it was shown that in cer- 
tain conditions this effect vanish. Besides this, for small 
aK, any such skew scattering or side-jump like current, 
even if present for non-uniform Rashba coupling, must 
again be quadratic in an, and is negligibly small. 

We now analyse the role of intrinsic spin-orbit interac- 
tions that contain only Gz Pauli operator for the spin cur- 
rent. The Hamiltonian for a single current-contributing 
conduction band reads 



"Hso == azUzA{k^,ky), 
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where A{kx:,ky) is odd with respect to the time reversal 
symmetry. In the case of interest for us, A{kx,ky) is 
defined by terms in Eqs.( 105|106|107 ). This Hamiltonian 
results in the spin-orbit term in electron energy: 
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We begin with a question whether the terms of this 
symmetry can contribute to the side-jump group cur- 
rents. These currents originate from two sourses: the 
"spin dipole" fi given by Eq.(59), and the phase of a 
scattering matrix element. It is easy to see that the 
spinors entering the electron wavefunctions and defining 
VL in this case are just [0, 1] or [1, 0], and are momentum- 
independent. Then the cooridnate contribution associ- 



ated with Vl in the side jump given by Eq.(58) vanishes. 
Furthermore, the "intrinsic" current related to curlVfi^ 
vanishes as well. Moreover, if a scattering potential it- 
self does not contain a spin-orbit term, the phase of the 
matrix element of scattering is spin-independent (and is 
due to different momenta in the initial and final states 
of a scattering process). Then, the phase part of the 



side-jump given by Eq.(58) also vanishes 



We next check if any contribution may appear as a 
result of renormalization of the velocity operator (and 
the current operator) containing the spin-orbit correction 
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Such spin-orbit correction is capable to produce the spin 
current, which is defined by terms odd in az in the ex- 
pression 



j± =-eX!v±/±W- 



(139) 



Here /±(k) is the nonequilibrium density matrix in the 
presence of externally applied electric field, given by 
Eq. ( 52 ) , which is linear in velocity. It is easy to see that 



when the expression under the sign of sum is odd in a^ , 
it is also odd in k, and the sum vanishes. Similarly, the 
contributions in which the spin-independent component 
of the velocity operator (1381 defines the current (1391, 



and the spin-dependence comes from energy (137), van- 
ish in summation over k. 

We analyse now the possibility of the skew scat terin g 
asymmetry caused by the spin-orbit interaction (136 1. 



First, due to momentum-independent and real spinors. 



the standard skew scattering probability given by Eq.(47) 
vanishes. Second, higher-order in scattering amplitude 
terms in the scattering probabilities also vanish. Indeed, 
an asymmetric distribution function in the presence of 
electric field is linear in the current (velocity) operator, 
and it contains odd in k term independent of a^ and 
even in k term linear in az- The current operator that 
is averaged over the distribution function has exactly the 
same property, while the spin-dependent part of energy 
is linear in az and odd in k. Then if an expression for the 
current is odd in a^, which is required for terms (1361 in 



order to get currents of opposite spin projections to be 
opposite in sign, it is also odd in k, so that summation 
over the directions of k gives zero spin current. 

Thus we see that nonequilibrium spin currents in the 
case of Hamiltonian ( 136 ) all vanish. However, there is 



quite an interesting twist here related to the spin cur- 
rents in equilibrium. This effect was pointed out by E.I. 
Rashba in connec tion with the spin-flipping Dresselhaus 
and Rashba term^i^U Por such contribution, the accom- 
panying spin relaxation of the z— component of spin is 
strong, and therefore the condition (74) cannot be sat- 
isfied, so that no spin-electric stripes occur. However, 
the "equilibrium spin current" caused by Hamiltonian 
(1361 is not accompanied by spin-relaxation. This ef- 



fect is given by the expression of Eq.( 139 ) with the func- 



tion /±(k) being the Fermi function. The linear in a^ 
spin current vanishes as an integral of the full differen- 
tial, but the cubic or higher power in a^, currents are 
present. These currents are very small, so that normally 
the residual spin relaxation, in particular due to nuclear 
spins, will preclude spin-electricity. Still there exists a 
theoretical possibility that in a setting with the spin- 
orbit Hamiltonian (136), parameter k ~ 0, and nuclei 
with spin 0, spin-electric stripes can arise. It would be 
curious to develop a material with these characteristics. 
However, we have no suitable candidate at present, and 
in what follows we will discuss only spin-electric stripes 
caused by the external electric field, i.e., the electric cur- 
rent flowing through the sample. 

As follows from the discussion above, the contribu- 
tions of all spin-orbit interactions except the principal 
term of Eq.(fl]) to non-equilibrium spin currents are ei- 
ther vanishingly small or vanish altogether. Therefore, 
our consideration of spin-relaxation as the only effect of 
the spin-flipping terms in a microscopic Hamiltonian on 
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the constituitive equations of the system Eqs. ( 62|63 ) is 
justified. 



C. Residual spin relaxation: electron-nuclear 
interactions 

In our quest to eliminate as much spin relaxation as 
possible we inevitably arrive to the point when the spin- 
orbit source of z— projection spin relaxation is reduced 
significantly, and even eliminated in an ideal structure 
with equal potential confinement of conduction and va- 
lence electrons in z-direction, with no change of mate- 
rials parameters accross the structure, while the spin- 
orbit terms giving electric field domains are intact. How- 
ever, as we are discussing III-V semiconductors as the 
most probable experimental setting, we must remember 
that these materials contain rather heavy nuclei in the 
host crystalline lattice, which are the source of hyperfine 
contact interactions of electron and nuclear spins. The 
contact interaction is usually the strongest electron spin- 
nuclear spin interactions, and potentially can result in 
electron spin relaxation. 

In the problem of spin Hall current, we are dealing with 
delocalized (free) electrons. Hyperfine coupling has never 
been considered as a viable spin relaxation mechanism 
for free charge carriers (as opposed to spin relaxation 
of localized electrons). For free electrons, the relaxation 
mechanisms associated with spin-orbit interactions are 
usually by far the most efficient, with relaxation rates 
5-6 orders of magnitude higher than for spin relaxation 
due to hyperfine coupling. However, because we aim at 
reducing (eliminating) spin-orbit interactions spin relax- 
ation mechanisms, we have to analyse the electron spin 
relaxation induced by hyperfine coupling. 

The hyperfine interaction of el ectr on and nuclear spins 
is described by the Hamiltonian,'i22l 
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where vq = eh/2mQC and vn = eh/2MpC are the electron 
and nuclear magnetons, mo and Mp being the free elec- 
tron and proton mass correspondingly, go ~ 2 is the free 
electron g-factor, index i enumerates the types of species 
of nuclei present in a crystal unit cell and n enumerates 
all nuclei of a given type in a sample, gi is the nuclear 
Lande-factor and li is the nuclear spin operator for the 
nuclei of i-type, i?i^„ is the nuclear coordinate. 

Spin relaxation of electrons due to hyperfine interac- 
tions comes from the following two mechanisms. One 
is the direct mutual electron-nuclear spin-flip scatter- 
ing. For the interaction given by Eq.(140), it is impor- 



tant that electron density is non-uniform over the crystal 
unit cell. For the conduction electron wavefunction in 
the ground state of an infinite quantum well Vl'fc(r_L) = 
(l/-\/^)w(r_L)expik_L • rj^(-\/2/(i) cosTTz/d, where u(rj^) 
is the Bloch amplitude periodic with the lattice constant. 



d is the width of the quantum well, and S is the 2D sam- 
ple area, we obtain scattering matrix element 
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where Ci^n = (2/d) cos^ {it Zi^n/d) is the coefficient arising 
from the quantum well confinement, Vc is the unit cell 
volume, Zi^n is the nuclear coordinate, and the hyperfine 
constant Ai is defined by 
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The spin- flip scattering probability is given by Eq. ( 96 ) 



with the substitution of the spin-orbit scattering matrix 
elements by Eq.(141|, with inclusion of averaging over 
initial states of the nuclei and summation over final nu- 
clear spin states. Averaging over coordinates Zi^n, gives 
the probability of the flip of the z— projection of spin of 
the quantum well electrons due to hyperfine interaction 
and the corresponding spin relaxation rate l/rj^fs- 
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where rii is the density of nuclei species of i-type. 

The second mechanism of electron spin relaxation 
comes from the effective magnetic field Bn acting on 
electrons spins due to hyperfine interaction ( 140 ) : 



Bn — 
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where angular brackets denote averaging over all nuclear 
wavefunctions and g* is the effective electron g-factor in 
a crystal. If this field were static, it would induce no 
electron spin relaxation. To analyze its dynamics, we 
first observe that there is also an effective magnetic field 
-Bg due to hyperfine interactions that acts on each nuclei 
of i-type due to electrons: 
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The effective field B^ ^ Bn, in particular because the 
nuclear magneton and spin splitting is 3 orders of magni- 
tude smaller than that for electron. We now observe that 
in the absence of the external magnetic field, and other 
sources of nuclear spin polarization, the effective field B^r 
averaged over times t > H/vNgiBl must be zero on sym- 
metry grounds. On a smaller time scale, nuclear spins 
retain some polarization caused by electron spins, and, 
in turn, the electron spin is subject to a frozen fluctua- 
tion of magnetic field arising from nuclear spins. These 
fluctuations of hyperfinc-induced effective magnetic field 
acting on electron spins are described by the dispersion 
of the nuclear field distribution, see, e.g., 1^221. Assuming 
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the nuclear spin directions independent, for nuclei acting 
on quantum well electrons, this dispersion is given by 
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Substituting the sum over n by an integration, we obtain 

(147) 
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where the summation runs over different types of nuclear 
species in one unit cell. This dispersion is very small 
for delocalized electrons. The corresponding time scale 
becomes of order of the characteristic tim es fo r nuclear 
dipole-dipole relaxation time, t^^ ~ io-4^i02| Nm^lear 
dipole-dipole interactions, in contrast to hyperfine inter- 
actions, do not conserve the total spin of electrons and 
nuclei. Then, the nuclear dipole-dipole interactions, as 
a result of flip-flop of the nuclear spins, randomise the 
effective magnetic fields acting on electronic spins due to 
hyperfine interactions, so that on the scale of T^d, the 
nuclear fields seen by electronic spins are rapidly chang- 
ing. If the characteristic precession frequency of electrons 
spins is fiBg*\/^B ^ ^■^dd^' ^^^ electron spin undergoes 
random rotations by an angle of characteristic magnitude 
/is5* i/A^Tdd/Ti. Such diffusive motion of the electron 
spin results in spin relaxation characterized by the relax- 
ation rate 
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The shorter Tdd, the faster the electron spin relaxation. 
Therefore in contrast to direct scatttering of electron 
spins by nuclear spins due to hyperfine interactions, i.e., 
coUisional broadening, the dipole-dipole nuclear interac- 
tions result in motional narrowing. In an experimental 
setting with almost completely suppressed spin-orbit in- 
duced spin relaxation, both collisional broadening and 
motional narrowing due to nuclear spins can be impor- 
tant factors in weak relaxation of z— projection of the 
electron spin. 



IV. SPIN-ELECTRIC STRIPES IN 110 GROWN 
QUANTUM WELLS WITH LIMITED SPIN FLIPS 

We now calculate the magnitude of the applied 
electric field needed for the observation of the spin- 
electric domains, and the magnitude of the transverse 
electric field in domains. We consider [110] grown 
Ino.52-AloAs^'\s/InP/Ino.52AlQAsAs heterostructure, and 
take K = 0.33^. This is the lowest k in known III-V sys- 
tems. The Kane parameter in InP is 1.36 x 10~^^erg-cm, 
the bandgap is Eg = 1.35eV, the effective mass is 
m — 0.07377i(p^and the spin-orbit splitting of the valence 
bands Fg and r7 ia A = O.lleT^. This gives spin-orbit 
constant a = 1.92 x 10~^^cm^. 



We condider a rectangular quantum well with 
300Awidth, in sample with the planar dimensions 94/i 
m xl cm. We take a strongly disordered quantum well 
with well doping by donors, and compensating doping 
by acceptors, at a total level of donor and acceptor im- 
purity density rit — 2 x 10^^ cm^^. In addition, two 
symmetrically positioned remote doping layers provide 
charge carriers, defining the total carrier density at a 
level n = 5xl0^^ cm~^ in all settings we are going to 
consider. When donor density n^ exceed acceptor den- 
sity Ha, donors provide additional electrons in the quan- 
tum well, and weaker remote doping is needed to get 
this level of total density. In the setting with acceptors 
more numerous than donors, remote doping layers need 
both to compensate the action of acceptors that take elec- 
trons from donors and to provide the total carrier density. 
This can be achieved by taking different separations of 
(5-doping layers (setbacks) from the centerplane z = of 
the quantum well. 

At large setbacks, the dominant scattering mechanism 
are quantum well impurities. We take scattering to be 
short-range but assume that it mimics shallow donors 
and acceptors capable of providing or accepting one elec- 
tron. Scattering amplitudes for such donor and accep- 
tor scattering have the same magnitudes and opposite 
signs. We consider three different settings corresponding 
to cases I, II and III of the Sec. IA: Case I corresponds 
to total compensation, n^ = 10^^ cm~^, n^ = lO^'^ 
cm~^, where n^ and Ua are the donor and acceptor den- 
sities correspondingly. In the case II, the donor doping 
level is ri(j = 1.1 X 10^^ cm^^, and the acceptor doping 
level is n^ = 0.9 x 10^^ cm^^. In the case III, quan- 
tum well acceptors ovcrcompcnsate quantum well donors: 
Ud = 0.9 X 10^^ cm^^, Ua = 1.1 X 10^^ cm^^. Remote 
doping provides charge carrier electrons in all three cases. 

For case I, we consider the stripes first neglecting the 
spin relaxation of the z— component of spin, and then see 
how spin relaxation modifies the picture. The 5-doping 
layers setbacks are / = 11 00 A, p roviding 7i = 5xl0^^ 
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cm~ electron charge density^^^, which results in Fermi 
energy ep = 16.5meV. The used setbacks suppress both 
scattering by smooth random potential and the contri- 
bution of long range-fluctuations into spin-orbit inter- 
actions. The impurity momentum relaxation time in 
this setting is t = 0.4 x 10"^"^, and the mobility is 
/i = lO'^crn^ /V ■ s. In this case of total compensation, 
the the spin current is deflned by side-jump-like contri- 
butions of Eq. ( 57 ) only. The parameter characterizing 
the strength of disorder in this setting is kpi = 2, and 
the transverse (spin Hall) mobility ryi = 0.6cm^/(l^ • s). 
A periphery stripes of 15/im width and transverse elec- 
tric field Ex = 2.4:V/cni, and a central stripe of the width 
25 — 64/iTO are induced in such low-mobility sample by 
an electric field Ey = AkV/cm. This gives energy ac- 
quired by an electron in electric field on the scale of the 
mean free path equal to eEi = A.2meV or ^ 0.25i?i? 
This would lead to certain heating, but nevertheless, the 
observation of such stripes in InAlAs/InP/InAlAs het- 
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erostructures is feasible. Designing heterostructures with 
lower K will allow to decrease electric fields required for 
the observation of the effect and thus to reduce heating 
substantially. The spin polarization, electric field and po- 
tential distribution neglecting spin relaxation are plotted 
in Fig. [5) Fig. [6] and Fig. [2] correspondingly. 

The choice of parameters made above is consistent with 
Eq.([74| and has a comfortable edge. Indeed, the calcula- 
tion of the spin relaxation rate for the z— component of 
spin shows that two mechanisms are dominant; First, 
weak residual spin-orbit interactions due to spin-orbit 



constants (112), (114) result in spin-orbit s catte ring off 
quantum well impurities described by Eq. (124), which 
leads to 1/tss = 5.47 x lO^'s^^. Second, the biggest spin 
z— component relaxation comes from direct scattering of 
electron spin by nuclear spins, due to hyperfine interac- 
tions. For InP, In nuclei has spin 9/2, gin — 1.11, and the 
corresponding density of the Bloch electron at the nuclei 
is|uc(0)|2 = 7.63x10^5 cm-' 105. For P nuclei spin is 1/2, 
gp = 2, and |ue(0]p = 3.26 x 10^^ cm^^ ^j^g relaxation 
rate given by Eq.|l43| leads to 1/tns = 1-297 x lO^g-^ 
The spin relaxation due to random Rasbha terms is ~ 200 
times smaller than that due to spin-orbit scattering, be- 
cause of the low mobility and large setbacks. The spin 
relaxation due to fluctuations of the hyperfine field in the 
presence of nuclear spin dipole-dipole scattering is very 
small due to massive sample. The total spin relaxation 
rate is l/r^ — 1.75 x lO^s^^ results in spin relaxation 
length Ls — 95.3/iTO. Thus, in the chosen structure the 
strip is narrower than Lg , which is more than 3 times big- 
ger than the half-width of the central stripe, and more 
than 6 times exceeds the width of the periphery stripe. 

The big spin relaxation length and long times are the 
consequences of weak spin-orbit constant in InP, struc- 
ture of offsets in the InAlAs/InP/InAlAs heterostruc- 
ture, "spin-orbit engineering", particularly, [110] crys- 
tallographic orientation of the quantum well, symmet- 
ric confinement, large setbacks for J-doping layers, and 
a wide quantum well. For comparison, the spin relax- 
ation rate due to spin-orbit scattering in the [110] GaAs 
quantum well, in an identical symmetrical structure, is 
IJTss = 0.65 X lO^s^i. In a lOOAwidth quantum well, 
I/tss = 0.6 X 10*5^1, and the corresponding spin relax- 
ation length is Ls = 4.7/ir7i. 

For case I, we calculate how the spin relaxation changes 
the picture of the periphery and central stripes. The half- 
width of the central stripe becomes bs — 39/im, and the 
width of the periphery stripe is now 8fim. The spin po- 
larization, electric field and potential distribution across 
the sample are shown in Fig.[T9J Fig. [20] and Fig. [21] cor- 
respondingy. A small electric field linear with coordinate 
is induced in the center reigion (Sec. IIC). This electric 
field, however, results in sizable build-up of the potential 
in the central region over the increased half-width of the 
central stripe. Still, placing the inside-sample contact 
for voltage probe within 8 microns from the edge, one 
would observe almost a linear increase of voltage with 
the increase of external electric field, characterizing the 




FIG. 19: Spatial dependence of spin polarization on the co- 
ordinate transverse to the flowing electric current in the pres- 
ence of spin relaxation of the z— component of spin, case I. 
All parameters are discussed in the text. 
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FIG. 20: Spatial distribution of transverse electric field in the 
presence of spin relaxation, case I. 

periphery stripe, as opposed to markedly quadratic de- 
pendence in the central stripe. 

We now discuss parameters of the system for cases 
II and III. Here the electron mobility is the same as 
in case I. Because the acceptors do not fully compen- 
sate donors, a skew scattering current is present. In 
the low density samples that we consider, the side jump 
and the skew scattering mobilities have equal nragni- 
tudes in the absence of compensation, n^ — 0. For 
Ind — nal/ind + TT-a) = 0-1, as in cases II and III, the skew 
scattering is just 10% of the side jump-like currents. We 
remark that the spin drag effect that decreases the skew 
scattering contribution*^, for small density of charge car- 
riers (4 times less than in''^) and temperature T=10 K , 
in our setting reduces the skew scattering by just ^ 2% 
of its value. The electric field and potential distribution 
neglecting spin relaxation are plotted Ey = AkV/cm for 
case II in Fig. J7] and Fig. J8] and for case III in Fig. ]9] and 
Fig. 10 correspondingly. The transverse electric field in 
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FIG. 21: Potential profile in the presence of spin relaxation, 
case I. 

the periphery domain is E^ — 2.16V/cm in case II, and 
Ex = 2.64V/cm in case III. In these cases, a linear in co- 
ordinate electric field in the central stripe arises already 
in the absence of spin relaxation. That means that at in 
the presently available setting it would be rather difficult 
to distinguish side-jump-like and skew scattering contri- 
butions by measuring the distribution of electric field and 
potential. However, if a system is engineered with bigger 
a and smaller k, or weaker constants of the hyperfine in- 
teractions, such an experiment would be feasible. We also 
note that case IV with dominant skew scattering or case 
V with skew scattering opposite in sign to side jump and 
exceeding it ^ 1.5 times in value requires the values of 
kpi ~ 10 that result in significant of eEy£ ~ Ep. These 
cases may become feasible in systems with stronger spin 
current and weaker spin relaxation. 

The electric field in both in stripes can be observed 
in experiments similar tcP^. A voltage between internal 
contacts in the center region and contacts on a sample 
periphery should be measured. A 100% spin polarization 
at the edges, opposite in two pperiphery stripes, can be 
measured in optical experiments similar tcP^. Furthemore, 



as Fig. 21 indicates, a build-up of electric potential and a 
significant spin polarization can be observed even if the 
periphery stripes are not yet developed, in a narrower 
sample. 



V. FUTURE WORK. CONCLUSION 

Considering the charge and spin currents, we restricted 
ourselves to currents linear in the gradient of density and 
electric fields. We considered zeroth approximation in 
charge-electric field distribution and calculated the first 
order modification of charge distribution, but did not 
consider the corresponding changes in accordance with 
the consituitive equations of the system. However, the 
approximations made capture the principal features of 
the spin-electric state. The current manusript is pur- 
posedly analytic, and all figures represent the analytic 
equations presented in the paper. Future work will ad- 
dress numerical analysis of time-dependent picture of the 
phenomena, self-consistent calculation of charge-electric 
field distribution together with the constitutive equa- 
tions, and hydrodynamic approach. Also, in this paper 



we restricted our consideration to conduction band elec- 
trons. The account of spin-electric stripes in systems 
with charge-carrier holes will be presented shortly. 

In conclusion, we the reiterate the principal findings of 
our work. When the electric current is flowing through a 
two-dimensional (2D) uniform conductor, the spin-orbit 
scattering results in the two spin-electric stripes near the 
two opposite boundaries of the conductor parallel to the 
direction of the current. In each periphery stripe there 
is an electric field transverse to the flowing current. The 
directions of the electric flelds in the two stripes are oppo- 
site to each other, and their magnitude is the same. If the 
relaxation of the spin component perpendicular to the 2D 
plane is negligible, the magnitude of the electric fleld is 
constant throughout the stripes. In the stripes, electrons 
are fully spin polarized, with opposite spin orientations 
in the two stripes. The electric flelds accompanied by 
a 100% spin polarizations in the periphery stripes are 
the consequences of the inability of the spin relaxation 
to limit the spin accumulation when a spin current flows. 
Spin polarization is limited only by its maximal (mini- 
mal) value. This leads to Anomalous Hall effect in the 
periphery stripes, so that charge carriers with the two 
opposite spins stream to the respective boundaries. It 
is then the Coulomb interactions that limit the streams 
and give rise to the electric flelds. 

The two periphery stripes are separated by the center 
region (the third stripe) of the 2D plane. On the center- 
line, the electron spin polarization is zero. The picture of 
the spin polarization in the central stripe is almost like 
in the conventional spin accumulation, except that when 
the spin polarization grows from the center-line to the 
boundary between the central and the periphery stripe 
(or falls for the opposite spin direction) in absence of 
spin relaxation, it reaches the maximal ±1 value at such 
boundaries and cannot grow any further into the periph- 
ery stripes. The magnitude and direction of the electric 
fleld in the central stripe in the case of negligible spin re- 
laxation depends on the relation between the microscopic 
mechanisms of the bulk spin current. Experimental mea- 
surement of the voltage caused by the electric flelds in 
the periphery stripes would constitute a direct electric 
measurement of the spin and AHE currents. In the case 
of reaching the limit of negligible spin relaxation, mea- 
suring the spatial distribution of the electric fields can 
provide a useful tool for experimental separation of the 
skew scattering and side-jump like effects. 

When a weak relaxation of the spin component normal 
to the 2D plane is present, the spin polarization and the 
electric fleld in the periphery stripes experience a small 
change from the periphery towards the center region. At 
strong spin relaxation stripes do not arise. The division 
of the sample into three stripes remains meaningful in 
the presence of weak spin relaxation: at the boundaries 
between the preriphery stripes and the central stripe the 
electric fleld is a steep function of the transverse coor- 
dinate. Also a 100% spin polarization is present at the 
edges of the sample. If a single periphery stripe becomes 



a part of a separate circuit with the electric current flow- 
ing, it is possible to transfer this spin polarization for the 
purpose of applications. 

Favorable experimental settings, in which an electron 
spin relaxation of the spin component transverse to the 
2D plane is suppressed but the spin current is not, are 
discussed. The role of the remote doping by donors and 
doping in the quantum well by donors and acceptors is 
studied. In a consideration of the spin-orbit effects for 
2D electrons, a 3D short-range and long-range random 
potentials are taken into account. The role of the bound- 
ary conditions and the role of the difference of parame- 
ters in the barriers and the quantum well are investi- 
gated. All the spin-orbit spin-flipping terms can be elim- 
inated from the 2D Hamiltonian by engineering a confine- 
ment potential for conduction and valence electrons. In 
the proposed experimental setting of electrons confined 
to InAlAs/InP/InAlAs symmetric double heterostruc- 
ture, it is possible to suppress the spin-orbit induced re- 
laxation of the normal component of spin to such an ex- 
tent that hyperfine interactions of spins of delocalized 
conduction electrons and nuclear spins play a significant 
role in electron spin relaxation. 
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Appendix A 

Here we present the details of a solution of the elec- 
trostatic problem, and find the electric field in the 
whole range of x and the distribution of charge den- 
sity across the 2D sample consistent with the solution 
of constituitive equations (pi) for the electric field in do- 
mains. Rather than deriving the function of complex 
variable E(t) of section LB from the complex potential 
for 2D electrostatic problem found by using the rele- 
vant Schwartz-Kristoffel transformation, we will derive 
the Kramers-Kronig like dispersion relations for some 
auxiliary analytic function, which is directly related to 
E{t). 

We consider a contour C in the complex plane t shown 
in Fig. [22j and a conformal mapping descibed by the 
function 
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The function f{t) has three extrema with the derivative 




-X, 



-bi -X, 



bi 



Xi 



FIG. 22: Contour in the complex plain t. The functions f{t) 
and 1/ fit), are analytic inside and on the contour formed by 
real axis, semicircular detours of infinitesimal radius around 
singular points ±a, ±bi, ±&, around points of extrema of fit) 
at ±^1, ±X2 and 0, and around infinite-radius semicircle. 

im = at 



h = (a' + ^{a^-bl){a^-b^))^/^ \t\ > a 

t2 = (a^ - v/(a2-62)(a2_fe2))i/2 fo < |i| < 5^ 

t = (150) 

Thus, the functions f{t) and l/f{t) are analytic in 
the upper half plane of the complex variable t and 
on the real axis (with the exception of points t = 
±ti, ±a,±6i, ±^2, ±^, 0). The complex function E{t) 
is analytic in the same region of the complex plane 
as well, and so is the auxiliary function Bit)/ fit). 
Furthermore, from boundary conditions E{t) given by 
Eqs.( |35|36|37|38| and from Eq.([l49| 



lim E{t)/f{t)^Q. 



1*1- 



(151) 



Then, from the integral Cauchy formula, analytic prop- 
erties of E{t)/ fit) and Eq.( |151 ), we have 



E{t)/f{t) = -P 



E{u)/f{u) 



t 



du. 



(152) 



where P indicates the Cauchy principal value. 

The real and imainary parts of E{x)/f{x) obey the 
dispersion relations 



1 



-P 



Re[Eit)lf{tl 
Im[E{t)/f{t)] =-Ip j 



Im[E{u)/f{u) 



— 00 

00 



u — t 

Re[E{u)/f{u 



t 



du (153) 
-di(.154) 



The function f{t) with the branching points ±a, ±61 and 
±6 is alternatively real and imaginary in the intervals 
|a;| > a, &i < |a;| < a, 6 < |x| < 61 and |a;| < 6, {x is the 
real value) : 



f{x) = sign{x) 



{x'^-bf){x^-b^) 



if \x\ > a 



fix ± zO) = T'' 



(x^-bj)(x^-b'^) 



if bi < |x| < a 



f{x) = s^gnix)^/^^^^ ii b < \x\ < b, 
/(x±zO) = ±z./i^gp^ if\x\<b (155) 
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Using dispersion relations (154), boundary conditions 
Eqs. ( 35|36|37 ), and Eq.( 155 1 with complex values taken 
in the upper half plane, we find unknown electric fields 
E{\x\ > a) and E{b < \x\ < bl) and charge densities 
n{bi < \x\ < a) and n(\x\ < b). For an unknown electric 
field, i.e., the electric field in the intervals of x, where it 
is not set by the boundary conditions, we have 



E{x, \x\ > a) — sign(x 

2 Eg fo- sds 



(2;2_fc2)(j,2_(,2) 

(a;2-a2) 



'4iVo/, 



bl sds 



(s2-6^)(s2_f,2) 

{bj-s^){s^-b^ 



(156) 



E{x,b < \x\ < 61) = sign{x) 

2E| 






2Eo p 

7T Jbi 



sds 



(s2-62)(^2_b2) 



'4^0 J {,2 s'^-t^y a2-62 



(b^-s2)(s2_j,2) 



(157) 



The charge density in the regions, where it is not de- 
fined by the constraint, is given by 



n{x,b, < \x\ <a) = -^ ' ^J>\,^ 

Eq r"- sds 

71-2 Jbi S^~t 



(s2-62)(g2_ft2) 



2^ r 

TT J02 



n{x,b < \x\ < bl) = 

-Eq f sds 

71-2 Jbi S^-t^ 



) No fbl sds 



{bl-s'^){s^-b^) 

?23^2 

(fc2-^2)(2;2-62) 
(a2— a;2) 

(s2-62-|(g2_fc2) 



(158) 



(62-s2)(.s2-62) 
^23^2 



(159) 



Using Eqs.(154| 156 1,(157 1, 



check that the electric field and the charge density in the 



1581 and (159), one can 



regions where they are set by Eqs.(|35[), (36) and (37) 



are given precisely by these values. Also, the complex 
function E{t) is indeed analytic and satisfies the condi- 
tion (38). Evaluation of integrals in (156), (157) and Eqs. 
(35|36), yields the resulting electric field profile given by 



Er. 



2-Eo 



^ i In I -h arctan ^^^ 



^/a' 



v^ 



2ttNo 



a^-b^ 



E„ 



signix) if \x\ > a 
Ex = EQsign(x) if hi < |a;| < a 
arctan ^^352 J 



(160) 
(161) 



2Eo 



W 



In 



a^-bl 



Va2-a;2+^a2-b2 

ln2J sign{x) if & < |a;| < 61 (162) 
Ex = if |x| < b, (163) 



where 



y 



^{bl-b^)/ia^-b^). 



(164) 



At large |a;| — >■ c», the asymptotic expansion of the 
electric field gives, on the one hand. 



|i?.l(W^oo) = ^?Hlf?^ln^+2.iVo 



(165) 
On the other hand, the magnitude of electric field at 
|a;| — )■ 00 is related the electric charge in the sample (see 
e.g.R, 



\Ex\{\x\^^)^'^=0, 



(166) 



where q is the linear charge density q, i.e., the charge per 
unit length in y— direction, which vanishes for an electri- 
cally neutral sample. 



From Eqs. (165 166), we determine 61 



bi=b 



8(a2 



s26 



exp' 



2n^No 

En 



(167) 



We see indeed that the constraint on charge density turns 
out to be important almost at a single point 61, which 
merges with b with an exponential accuracy. However, 
neglecting the constraint would result in logarithmic di- 
vergency of the charge density. Evaluating the integrals 



gives the charge c 



in Eqs. ( 158 ) and ( 159 ) taking the constraint into account 



ensity profile: 






^a2-bfK^2_fc2)_ 

(a2-3;2)(a;2-62) i" ^ 



c2_b2 



In 



^-b'i (Va2-a;2 + ^a2-62)2 



a2-62 



In 



V'^^^yb^b^ 



(Va;2-fc2+yj.2_f,2)ya2-b2 

if bl < X < a 



'No if b2 < X <bi 



No 



a^-b'i 



(168) 



(169) 



^« ' i In 



bj-x' 



(a^-bl)(b^-x^) , 
(b2-a;2)(a2_a;2) i^i ^ 



+ 



1 [W- 



In 



(Va^-a;2 + ya2-b2)2 



+ 



a^-b^ 



In 



^/^^^^/bf^ 



(Vb2-x2 + yj,2_f,2)ya2-b2 

if X < &2 



-No 



a^-b^ 



(170) 



Electric field and charge density profile are presented in 



Fig. 14 and Fig. 15 for case I parameters given in Sec. VI. 
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Appendix B 

Spin-orbit constants and interactions in various 

systems 

In this appendix we give the values of spin-orbit scat- 
tering constants. We will present a discussion of spin- 
orbit effects in metals, bulk semiconductors and 2D elec- 
tron gas in semiconductor quantum wells. Spin-orbit in- 
teractions have been treated in numerous monographs 
and textbooks. Nevertheless, the phenomenon we deal 
with in the present paper requires a special vantage point: 
how to suppress the spin-orbit effects resulting in spin re- 
laxation of the z— component of spin while keeping the 
part of spin-orbit interactions responsible for the spin 
Hall current intact. This is our primary focus here, which 
requires development of theory of spin-orbit scattering in 
quantum wells, as well as re-consideration and clarifica- 
tion of several facts about spin-orbit interactions, includ- 
ing intrinsic spin-orbit effects, e.g., of the Rashba type. 



Magnitudes of constants describing electron spin-orbit 
interactions in bulk systems 

Spm-orbit constants for relativistic electrons and electrons in 

metals 



In relativistic quantum mechanics, the spin-orbit scat- 
tering potential in the Schroedinger-Pauli equation arises 
as a result of spin-dependent admixture of positron states 
to electron states by coupling proportional to the speed 
of light and the electron momentum in the Dirac Hamil- 
tonian. The standard hand-waving explanation of the 
spin-orbit coupling is that a moving electron experiences 
the magnetic field vWU/c that couples to its magnetic 
moment. This explanation results in almost the same 
magnitude of the relativistic spin-orbit constant as the 
rigorous derivation from the Dirac equation, except it is 
a factor of two higheplsni. 

For ordinary metals with a broad conduction band 
and the Fermi energy ~ lOel^, the resulting spin-orbit 
constant for electrons has essentially a relativistic value 
a = h/4m?c^. Electron-electron interactions in a Fermi 
liquid may somewhat renormalize this constant, in much 
the same way they renormalize the effective mass of elec- 
trons in metals. Modification of the magnitude of the 
spin-orbit constant due to an admixture of states with 
different symmetry to electronic states (similar to the ef- 
fect that we will discuss below for semiconductors) is very 
small for metals. Most importantly, in simple metals, like 
Al, Cu, or Au, the spin-orbit interaction is isotropic, i.e., 
the relaxation of spins a^, Uy and cr^ is defined by the 
same time constant as in Eq.(97). In order to be able 



to suppress spin relaxation for one component of spin 
but to keep the terms responsible for the spin current, 
significant anisotropy of spin-orbit interactions is necces- 
sary, which is not the case for simple metals. We need 
the spin-orbit interaction given by Eq.(IT]) be dominant. 



while similar terms containing ax and Gy be small. This 
might be possible in a two-dimensional metal, but such 
structures similar to the 2D semiconductor systems have 
not been realized. The apparent reason for difficulties 
in designing metallic quantum wells, in which only one 
ground level of size quantization is occupied, is a very 
high electron density in metals. 



Spin-orbit interactions in semiconductors, including 
heterogeneous media 

In semiconductors, the spin-orbit constants are signif- 
icantly different from those in metals. The solution of 
the Schroedinger-Pauli equation in a periodic potential 
V{y) results in a group of bands, which are energetically 
close to each other. In particular, in the Kane modeP^ 
describing III-V systems, the conduction band \c) (group 
theory notation Fg) is separated from the top of the va- 
lence band iFg) by Eg just ~leV. The band Fg with the 
total angular momentum 3/2 is separated from electron 
states in the split-off valence band jFy) with total angular 
momentum 1/2 by the spin-orbit splitting A 

^ = -^^(-l(Vnr)xp),|z), (171) 

where functions (x|, (y|, {z\ are orthonormalizcd func- 
tions that behave like vector components x, y, z upon 
the symmetry transformations of the crystalline point 
group. The spin-orbit splitting in III-V systems can be 
smaller or bigger than the band gap. Furthermore, in 
much the same way as the Dirac electrons are coupled to 
the Dirac positrons by spin-dependent coupling, the con- 
duction band electrons are coupled to the valence band 
electrons, so that the Kane Hamiltonian T-L^ matrix ele- 
ments are given by: 



'^c,l/2;r8,l/2 ~ V ^^""^ 



K 



1 



fc_P 



■T/" 

'T-c,l/2;r8,-l/2 — /g 

^ - k 



''-c,l/2;r8,3/2 



''•c,l/2;rs,-3/2 

-w^ - -J: 

"c -l/2;r8,l/2 - /g 





k+P 

k,P 



njK 

'T-c.-l/2;r8 -1/2 



^c,-l/2;rs,3/2 — 
'^c,-l/2;r8,-3/2 — ^-P 



''-C,l/2;r7,l 



/2 



V3 



KP 



'^c,l/2;r7.-l/2 — \ o^^P 



'^c,-l/2;r7,l/2 — ^\/ 3^-^ 
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n 



K 

c,-l/2;r7,-l/2 






(172) 



where the Kane matrix element P — —{S\px\X), \S) is 
the orthonormahzed function that behaves hke a spher- 
ical function upon the symmetry transformations of the 
crystalline point group, \X) is the orthonormahzed wave- 
function that behaves like an x-component of the coor- 
dinate upon these symmetry transformations, and k± — 
(kx ± iky)/^/2. When the Kane parameter changes with 
the coordinate across the heterostructure, the order of 
operators kz and P{z) in Eq.(172| is important. Below 



we will discuss the choice made for the order of these op- 
erators in Eq. ( 172 ) and its consequences for the form 



of the conduction band Hamiltonian, the effective mass, 
the spin-orbit constant and for the boundary conditions 
at heterostructure interfaces. At this stage, it is impor- 
tant that such choice results in the corresponding form 
of the Hermitian conjugated terms of the Kane Hamilto- 
nian, e.g.: 



njK 

''•r8,l/2;c,l/2 



IpK 



(173) 



The spin-orbit interaction in the conduction band of bulk 
semiconductors arises in the third order perturbation the- 
ory. We shall see here that the magnitude of spin-orbit 
interactions is affected by a spatial dependence of the 
electron mass, P{z) and A(z). Then, both the second 
and third order perturbation theory terms are important 
for us. In the second order, the kinetic energy operator 
associated with the z— direction motion becomes: 



H 



kin 



ld_ 
3 9^ 



P{z) 



Eg{z) Eg{z) + Aiz) 



Piz] 



1 i„, ,_a 

dz 

(174) 

We see that our choice of ordering of the operators in 
the Kane model leads to a form of kinetic energy that is 
described exclusively in terms of the position-dependent 
Kane effective mass given by 



TO 

TOc(z) 



2to „, , 

-^P{z) 



2 1 

Eg{^)^ Eg{z)+A{z) 



Piz). (175) 



The total kinetic energy, assuming that the effective elec- 
tron mass changes along the growth drection only, is 

A choice of any other ways of ordering of P{z) and g| in 
the Kane operator would result in description of motion 
of electrons in z— direction by three varying parameters 
(P(z), Eg(z), and A(z)), instead of a single parameter 
mc{z). It is noteworthy that electrons in the conduc- 
tion band with varying effective mass have been often de- 
scribed in the literature by the Hamiltonian, see, e.g.p^^ 



n: 



kin 



d 



d 



dz (toc(z))i+2" dz 



(177) 



Therefore we see that the Kane model leads to such de- 
scription only in the case a = 0, which corresponds to 
our choice of ordering operators. We now consider the 
spin-dependent part of the electron Hamiltonian. The 
admixture of the valence band to conduction band elec- 
trons in the Kane model gives 



1 



[V, 



P^z) 



E-U,{r)-Aiz) 



X k] • (T, 



(178) 



where E is the electron energy, the first term comes from 
Fg bands and the second term originates from the split- 
off band. The spin-orbit coupling ( 178 ) is defined by such 



non-relativistic parameters such as the Kane matrix ele- 
ment P and the band gap Eg . It is related to relativistic 
spin-orbit coupling via spin-orbit splitting of the valence 
band A and vanishes in the limit of zero A. The mag- 
nitude of spin-orbit interaction critically depends on the 



ratio of Eg and A. An important feature of (|178|) is 



its dependence on the potential C/t,(r) that differs from 
potential acting on c ond uction electrons. This feature 
attracted attention irP^ and was discussed irP3 in con- 
nection with the consideration of the Rashba term in the 
2D electron gas. As we shall see, the story has much 
broader implications, and this feature is crucial for under- 
standing the constant of spin-orbit scattering in quantum 
wells, spin-orbit coupling due to long-range fluctuations, 
and the possibility to suppress both spin-orbit scattering 
and various Rashba effects. 



Effective 3D Hamiltonian for abrupt hetero junctions. 
Boundary conditions. 

Our primary interest are spin-orbit effects in symmet- 
ric quantum well structures with abrupt heteroj unctions. 
We will first obtain "3D" Hamiltonian describing such 
systems, and then derive the 2D effective Hamiltonian 
by averaging over the wavefunction in the ground quan- 
tum well subband. The potential affecting conduction 
electrons is given by Eq.(108l, the potential for Eg va- 
lence electrons is given by Eq.(108), and for r7 electrons 
by Eq.(lll). Assuming the coordinate-dependent parts 



of the potential affecting valence electrons to be much 



smaller than E"^ and E"^ 



Aon, we find 



H.o = -MVr(p2(^)t/.(r)) X k] • a[^ 



(B„)2 



(£™+A„ 



-iV,(p2(z))[^ 



P\z)5{A{z)) 

1 A 



[kxa]^. (179) 



We observe that the spin-orbit Hamiltonian includes 



three contributions. The first term in ( 179 1 is a gen- 



eralization of standard 3D spin-orbit interactions 

H'° = aa- [VrC/„(r) x k]. 



(180) 
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with the constant a given by 



by 



a 



PV.l 



3 VK2 {E, + Ar 



(181) 



to a heterogeneous media with a varying Kane constant. 
We will use these equations to describe spin-orbit effects 
in the approximation neglecting all differences in mate- 
rials constants in the well and the barriers, except for 
Eg. Then the materials parameters in Eqs.( 180] 181) 
will correspond to their values in the quantum well. 

The spin-orbit coupling associated with the first term 
in Eq. (179) arises from coupling of the conduction and 



the valence bands in a manner similar to the appear- 
ance of the relativistic spin-orbit coupling from electron- 
positron coupling in the Dirac Hamiltonian. This is a 
third order perturbation theory term: the first stage of 
the process comes from coupling of the conduction and 
the valence bands given by Eq.(172), the second stage is 



an action of the varying potential Uv{r) (including off- 
sets) acting on the valence electrons, and the final stage 
couples valence bands back to the conduction band. The 
specific form of the first term (179), i.e., the order of 



differential operators and P{z) is defined by our choice 
of ordering of these operators in the Kane model. The 
second term in ( 179 ) comes from the difference of the en- 
ergy separation between valence bands Fg and Fy in the 
quantum well and in the barriers. Finally, the third term 
arises in the second order perturbation theory and is as- 
sociated with the difference of the Kane matrix element 
P in the barriers and the well. The specific form of this 
term is also the consequence of our choice of ordering of 
the operators in the Kane model. 

We now discuss how the ordering of operators in the 
Kane Hamiltonian affects the boundary conditions for 
the wavefunction and spin-orbit interactions. As dis- 
cussed i]J25|li09|[Tiol^ the general way to order operators 
in the Kane model is to write down the matrix elements 
in the form, 






/2-XsM2 



P^KP' 



(182) 



with < a < 1- The multiband wavefunc tion in the 
Kane model is constructed, according to SuriaSH^ in the 
form 



V' = u|5)+v|R), 



with 



and 



"1/2 
U-l/2 



Vl/2 
V-1/2 



(183) 



(184) 



(185) 



where u; and vn {I = ±1/2, i = x,y,z), Rj, = X, R.y — 
Y , and Rz = Z. Then the continuity equation is given 



— (luP -I- Ivp) -I- Pdiv(u+w + UV+) = 0. (186) 

ot 



The conserved flux is 



I = P{u+v 



(187) 



and P^u and P^ 



are continuous at the boundaries. 



where z is the normal to the boundaries. Then, relating 
u and V from the Kane equation T-L^'i}} = Eip , and using 
the continuity of P^~°'Vz, we obtain continuity of the 
quantity 



p2 



E+E„ 



E+E„ 



du 



A[a-k].M 



{E+Eg){E+Eg+A) 



(188) 



at the boundaries. At a ^ 0, u is not continuous. 
However, there is no contradiction here with general 
quantum-mechanical principles, because, at an abrupt 
heterojunction, the total wavefunction constitues a mix- 
ture of the wavefunction of many bands, and only a total 
wavefunction is subject to continuity and flux continuity, 
rather than its individual components, such as u. An ex- 
ample of boundary condition in the Kane model, in w hich 
u is continuous but {du/dz)/mc is not, was devired irP^. 
Remarkably, our choice of ordering in the Kane Hamil- 
tonian with a = at _E = 0, as follows from Eq.(188l, 



results in the continuity of u and {du/dz)/mi, themselves, 
i.e., in the natural and traditional boundary conditions 
in the quantum well. We also note that the model consid- 
ered here for coupling of conduction and valence electrons 
is a "pure" Kane model (a model with no bare diagonal 
terms quadratic in k with different masses in the conduc- 
tion and valence bands). Therefore the choice a = is 
fully consistent, and there is no reason for the asymmetric 
choice a = 1, such as made in modified Kane modeil^^. 

We also note that for every choice of a, our principal 
conclusions about spin-orbit coupling will hold. We see 
that in order to reduce spin-flipping terms, materials in 
the quantum well and the barrier should have close values 
of materials constants. 

Concluding this part we remark that in semiconduc- 
tor physics, the flve band and seven band models are 
sometimes used, which are more elaborate than the Kane 
model. In these cases, it is possible to obtain contribu- 
tions to spin-orbit coupling in which a and VrC/«(r) are 
coupled to an expression cubic (or higher power) in com- 
ponents of the momentum operator. These terms are 



much smaller than spin-flipping terms in (180), and alter 



the symmetry of the Hamiltonian in a minimal way, not 
essential for our purposes. 



Spin-orbit effects in the ground conduction band subband of 
2D electrons in quantum wells 



Having derived the 3D Hamiltonian ( 178|179 ), we now 
proceed to describing the spin-orbit coupling for 2D elec- 
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trons. We will temporarily ignore the effect of different 
masses, P and A in the barrier and the quantum well, 
which will be taken into account in Appendix C, and as- 
sume that the only feature of a heterostructure is the 
potential confinement characterized by the offsets in the 
conduction and valence bands. 



The effective Hamiltonian for quantum wells 

The most general potential describing charge carriers 
in quantum wells is given by 



where 



U,{v)^U{z) + V{x,v,z), 



U{z) = Usymiz) + Uasymiz) 



(189) 



(190) 



describes the potential confinement of charge carriers to 
two dimensions, and is uniform in a;, y. The asymmet- 
ric potential Uasym{z) can include an external electric 
field applied along z. The potential (|190|) must result 



in a confined state or states, so that no free motion or 
escape into the delocalized state along the z— direction 
takes place. The separation of U{z) onto symmetric and 
antisymmetric components is convenient in the case of 
the double heterostructure rectangular quantum well, in 
which components of U{z) are symmetric or antisym- 
metric with respect to the centerline of the QW at z = 0. 
The second term in Eq.(|48]) can describe charge carrier 
scattering, fiuctuating electric fields due to large-scale 
fluctuations in impurity density distribution, or the ef- 
fects of fluctuations of the quantum well width by one 
or two monolayers. In two dimensions, we assume that 
V{x, y, z) does not result in confinement in transverse {x 
or y) directions, so that the in-plane motion of charge 
carriers is almost free, and is subject to effects of scat- 
tering by V{x,y,z). 

The closest situation to an ideal two dimensional case 
occurs when electrons occupy just one level of size quanti- 
zation along the z direction. Actual transitions of charge 
carriers to higher size quantization levels at low den- 
sities, such that the Fermi energy lies below the bot- 
tom of the band originating from the first excited level 
of the quantum well, require energy conservation, i.e., 
phonons, and are negligible at low temperatures. Thus, 
scattering leaves electrons on the same ground size quan- 
tization level. However, as we shall see shortly, due to 
a dependence of scattering potential on z, an admix- 
ture of excited electron states in the confining potential 
Usym{z) + Uasymiz) to the ground state plays an impor- 
tant role. 

The ideal 2D spin-orbit interactions described by 
Eq.Q correspond to the case oiV{x,y,z) = V{x,y), i.e., 
when no dependence of random potential on the third di- 
mension takes place. As we neglect for now the difference 
of masses, P and A in barriers and the quantum well, our 
starting point here is Eq. (180). 



In general, the Hamiltonian ( 180 ) is not reduced to 
the Hamiltonian (fl]) in two dimensions. Each impurity 
has a certain location in a quantum well along z. From 
a symmetry standpoint, this provides a vector lowering 
the symmetry of the system and results in symmetry al- 
lowed spin-dependent invariants in the Hamiltonian of 
the system. In a rectangular quantum well of a dou- 
ble heterostructure, the potential of impurities is neither 
symmetric nor antisymmetric with respect to z = 0. As 
a result, the Hamiltonian responsible for electron scat- 
tering in general includes not only a^ but also a^ and ay 
and can lead to relaxation of the z-component of spin. 

We take a realistic potential of Eq.( 189 ) and go beyond 



traditional approximations of a scalar potential being a 
sum of two functions, one depending on the growth axis 
coordinate and the other depending only on the in-plane 
coordinate. Furthemore, scattering potential in general 
is not even a product of such two functions. We begin 
with the derivation of the spin-independent scattering in 
two dimensions. The Schroedinger equation is 

2 2 

[f^ + 1^ + Uix, y, z)]^{x, y, z) = E^ix, y, z). (191) 
2m Zm 

We represent the wavefunction in a form 



^(a;,y,z) = *^j,(z)$(a;,y). 



(192) 



where 



[^ + Usyra{z) + Uasym{z)+V{x, y, z)]*. 



,iz) 



Exy^>Xy{z) 



(193) 

is the Schroedinger equation for z-dircction at the fixed 
values of in-plane coordinates treated as parameters, 
E{x, y) is an eigenenergy, and ^2,^ (z) is the wavefunction, 

with the normalization condition J_ dz\'^xyiz)\'^ = 1. 
The 2D problem is then separated into the two parts. 
First we need to find E{x,y) and_^^^(z). Our inter- 
est is the ground state solution of (193), so E{x,y) and 
"^xyiz) from now on refer to the ground level. We in- 
sert (192) into (191), multiply its both parts by *Jy(z), 
and integrate over the coordinate z, thus obtaining the 
Schroedinger equation for 2D electrons 



2 

^2to 



E„j + P{x, y) + L^y\^{x, y) 



E^{x,y), (194) 



where P{x,y) = {p± ■ A^ -I- A^ • p±)/2m describes an 
in-plane motion in the presence of an effective vector po- 



tential Aj^ 



This vector 



iJZodzKyiz)d/dr^^,y{z). 
potential and the associated magnetic field B" = curlA"' 
vanish if the wavefunction "^^yiz) is real. For localized 
states in the absence of magnetic field and absence of 
additional non-trivial dependence of V{x,y,z) on spin 
(coming, e.g., from magnetic textures or magnetic impu- 
rities), this is indeed the case. Furthermore, in Eq.(194), 
Lxy = -2^/_!^dz**y(z)a2/ar2_*^j^(z) is an effective 
scalar potential. Both P{x,y) nd L{x,y) are induced by 
scattering modulation [Eq.(193)] of the wavefunction of 
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the confined state. These terms are small only if the mod- 
ulation can be considered slow in adiabatic approxima- 
tion. However here such adiabatic approximation is not 
essential. We note that it is easy to generalize Eq.(194) 



to an effective matrix Hamiltonian that describes several 
excited states I of the Hamiltonian ( 193 ) in addition to 



the ground state, and takes into account both diagonal 
and off-diagonal terms of the 2D Hamiltonian in the basis 
of wavefunctions '^l. {z). 



Spin-orbit interactions in the quantum well. Offsets 
contribution 



The effect of spin-orbit interaction on the in-plane mo- 
tion arises due to an additional term, S{x,y), in the the 
Hamiltonian Eq.(194) when the 3D spin-orbit term, e.g., 



Eq.(180), is averaged over the z-direction. 



S{x,y) 



dz^:y{z)H^°^,y{, 



(195) 



If 7J*° is strong by some reason, the spin-orbit term 
must be treated not as a perturbation (1951, but rather 
included into V{x,y,z). Then, L{x,y), P{x,y), and 
E{x, y) will include the associated spin-orbit effects. Here 
we use the pe rturb ative scheme that includes spin-orbit 
effects via Eq.(195). 



We now consider the spin-orbit interactions in a partic- 
ular case of scattering in a symmetric rectangular quan- 
tum well of a finite depth. This results in a single source 
of spin-orbit effects, the scattering potential. We will 
easily generalize our result to the case of an arbitrary po- 
tential confinement. Because we are interested in a scat- 
tering problem, in order to average H"" over the ground 
state of the Hamiltonian (193), we will use the pertur- 
bation theory. In linear in V{x, y, z) approximation, it is 
important to include two contributions: in (i) V{x, y, z) 
enters Hgo', in (ii) V{x,y,z) enters '^^yiz), while Hso is 



defined by Eq.( 190| or just by symmetric potential con- 
finement. In the first order in scattering 



^,yiz)=^^{z)+Y^ 



Von{x,y) 

Eq ~ En 



*"W, 



(196) 



where ^"(z) is the solution to the Schroedinger equation 
describing the z— direction only 



H^*„(z) 



^2m 



+ f/(z)]*„(z)-^„*„(z), (197) 



where index n includes both integers n = 0, 1, 2... corre- 
sponding to confined states of a discrete spectrum and 
a continuum of wavevectors corresponding to the con- 
tinuous spectrum. The summation over n in Eq.(196) 



includes integral over continuum of states and sum over 
the discrete spectrum. In the case when U(x,y,z) = 
U{z) + Vscat{x, y, z) is the sum of the quantum well con- 
finement potential and a scattering potential, for elec- 
trons with energy Ei on the ground level of the size 



quantization, the effective 2D potential E^y is 



E,+ dz'^*{z)Vscat{x,y,z)^^y{z). (198) 



-E.T 



Averaging over '^xy{z) here includes sum over all virtual 
excited states of the Hamiltonian (197), in accord with 
Eq.(p6l). 



Then the 2D spin-orbit term reads 

'^lo = aJZ,dz^lJz)[a,[V,U4r) x k], + 
V,[/,(r)[k X a], + [ax Vrt/„(r)],fc,]]*,,j,(z),(199) 

The first term in brackets under the sign of integral 
in ( |199| ) results in the spin-orbit interaction of the same 
symmetry as the ideal 2D case spin-conserving interac- 
tion ([I]). In the first order in Vscat this term is given 
by 



n 



2D 



a I dz%{z)[a,[VrVscatix,y,z) x k],'^a{z) 

(200) 



The two other terms in ( 199 ) potentially can lead to 



the relaxation of z— projection of spin. The second terms 
in brackets, in the case when a potential confining charge 
carriers to the quantum well is asymmetric, is known to 
result in the Rashba 2D term associated with this asym- 
metry. For Uv{r) due to scattering potential, this term 
produces spin-flip scattering, and for Uy{r) associated 
with fluctuations in S doping, a random Rashba term 
arises. Furthermore, the third term in brackets, which 
makes the Hamiltonian Hermitian also contributes to the 
relaxation of the z— component of spin due to scattering 
or the random Rashba interaction. 

We now analyze the second term in Eq.(199) con- 
tainiing VzC/„(r). Inserting Eq.(108) into Eq.(48), and 
Eqs.( |48|196|109[ ) into Eq.([l99|), we obtain the fol owing 
contribution to the 2D spin-orbit interaction: 



nl"^ ^ a(^J^^dz[%{z)\/,Vscat{x,y,z)^o{z) 



E 



K.o(a:,y) ,T,* 



n Eo-E„ 



**(z)V,[/„(z)*o(^)- 



+T{x,y) (201) 

Here the first term in the integral originates from the gra- 
dient of the scattering potential, and the other two terms 
originate from the gradient of the potential confining the 
valence electrons (even if this potential is symmetric), 
but i nclude the scattering potential via the wavefunction 
(196). In these terms, kj^ is an operator acting on wave- 



functions '^{x,y) in Eqs.(192 194) The last contribution 
T{x, y) reads: 



r(a:,2/) = a/_°°^dzE4([kxa] 
*;(z)V,C/„(z)vI'o(z) + vI-*(z)( 
V^t/„(z)*„(z) 



Vno{x,y) 
' Eo-E„ 



(202) 
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In T{x,y), the differential operators k_L in [k x a]z act 
only on Vno{x, y), which originates from the effect of scat- 
tering on the wavefunction of confined electrons. Cor- 
respondingly, due to T{x,y), the wavefunction $(x,y) 
of the transverse motion is affected by the operators 
idVnoix,y)/dx and idVno{x,y)/dy. 



We now transform the first three terms of Eq.(201| 



From Eq.(193), we have an identity 

V,C/(z) = -^[H^fc,], (203) 

where [A, B] ^ AB - BA. Therefore, from Eq.( |108|109 l 



v,c/„(z) = ^[H^fc,] 



(204) 



Then, e.g., the matrix element (V2C/(z)) , in Eq. (201) 
is given by 

(V.C/.(.))„„ = g(fc.)„„(i?o-i?„). (205) 

Inserting Eq.( |205| into Eq.( |201[ ), and using the identity 



O'^l'^^/nO V^^IOn^nO 



Kfr*l = (v^y'"^"*) 



(206) 



we arrive to the result 

lie 



nj2D 



a 1 



V,F^ 



'00 



[kxa],+T(a;,y). (207) 



Considering now the T{x, y) term, and applying the 
identitities of Eq. ( 204 ) and the identity 



E[(^^)o, 

n 

we obtain 



"^"0 1 _ (I, 



d 



/OO' 



(208) 



T(.,,) = «|((fc.|:^-'^*)oo-. - (^^""'"^-"^ 



'00 



(209) 
Finally, combining Eqs.(207 209) with the term in 



Eq.(199) containing fc^, for the spin-flipping terms of the 



2D spin-orbit interactions we obtain 



nj2D 



all 



Uv_ 
tJr 



[(v,y 



^^")ookx-]^ 



([axV^F-'^*].fc,)J. (210) 

These are spin-orbit interaction terms that lead to relax- 
ation of the z— projection of spin, due to asymmetric scat- 
tering potential V{x,y,z) in the rectangular quantum 
well due to different offsets in the conduction and valence 
bands. More generally, for the arbitrary shape of the elec- 
tron quantum confinement, the corresponding spin-orbit 
terms are defined by Eqs.( |201|202| ). From Eq.(p03|), it 
follows that such scattering-induced terms are non-zero 
only if the potentials acting on the conduction and va- 
lence electrons differ from each other. This generalizes 
a similar result for average Rashba spin-orbit interaction 
due to one-dimensional asymmetric potential^ to spin- 
orbit interactions due to the 3D scattering potential or 
due to long-range fluctuations. 



Appendix C. Effect of the difference of parameters 

in barriers and the quantum well on spin-orbit 

interactions in 2D systems 

We now address the contribution to spin-orbit effects 
associated with the change in the effective masses, Kane 
matrix element P, and the spin-orbit splitting A at the 
heterostructure interfaces. We will begin with the ef- 
fect of the kinetic energy term in Eq.(174). The pres- 



ence of this term leads to two contributions. First, ob- 
viously a decay of wavefunctions in barriers and oscilla- 
tory behaviour of the wavefunctions inside the quantum 
well are described by different masses, and this changes 
the magnitudes of matrix elements in Eq.(210). Second, 



an additional non-trivial contribution to spin-orbit terms 
arises due to modification of Eqs.( 203|204 ). Indeed, in 
the presence of mass difference described by the Hamil- 
tonian Eqs.( 174|175 1, the gradient V zU{r) is given by 



V,C/(z) = -^[H^fc.] + 



d 1 dm{z) 1 d 
dzm{z) dz m{z) dz' 



(211) 



and we now take into account the second ternr in this 



identity. Using Eq.(108 1091 we have 



V,C/„(r)=V,C/(r)^ 



(212) 



Then from Eq.(180) the 2D spin-orbit interaction due to 



mass difference in the quantum well and barriers is given 
by 



U, 



SV-so ~ cth {niB — mw) — 

-"^zKyidmk X a],-V,*,^(d/2) - 
Lm " m 

-V.*:^(-rf/2)[k X al-Vz^,yi~d/2)\ . (213) 



We note that the flux —'^z'^xyi^) is continuous at the 
interfaces for chosen boundary conditions. Importantly, 
as a result of scattering or due to the asymmetric po- 
tential, the flux values are different at z = — d/2 and at 
z = d/2. We also note that in the presence of scattering 
potential V{x,y,z), the operator [k x a]z that includes 
spatial derivatives over in-plane coordinates acts both on 
^{x,y) and on \E'a;y(±(i/2). The latter i n the first order 
in scattering potential is defined by Eq.(196|. The sum 



of the corresponding terms gives a Hermitean spin-orbit 
Hamiltonian. 

When only the ID asymmetric potential Uasym{z) is 
present, the spin-orbit interaction due to mass difference 
in the well and barriers contributes to the Rashba term: 

5Hl^ ^ a{mB - mw)%'y 



[^^,^^*(d/2)^Vz^{d/2)- 
iV.**(-d/2)iV.*(-d/2)][kxa], 



(214) 
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In a weak asymmetric potential, e.g., due to gate voltage and the constant 7 is given by 
applied to a symmetric quantum well or due to asym- 
metric doping of symmetric heterostructure, in the first 
order in the asymmetric potential we have 



^ 1/J 

^ 3\E^ E^ 



(218) 



^{z)^¥\z)+J2 



Eq — -E„ 



*"(Z), 



(215) 



A difference of the probabilities at the left and right het- 
erointerfaces, or a difference of products of the wavefunc- 
tion and its derivative over in-plane coordinate at the 



and inserting the second term into Eq.(|214|) results in a left and right heterointerfaces entering Eq.(pT6l) is due to 



Rashba term associated with potential asymmetry. 



For strongly asymmetric quantum wells Eq.(214| will 



contain the difference of fluxes at the interfaces of the 
effective mass change z = —d/2 and z = d/2. We note 
that the contribution to spin-orbit terms due to differ- 
ent masses will also take place in inversion layers. Such 
contribution can be calculated using general expressions 
for the spin-orbit terms Eq.(199) and Eq.(174). Strong 
asymmetry in inversion layers will lead to sizable Rashba 
term as compared to spin-orbit terms that conserve the 
z— projection of spin. 

We now consider the effects described by the second 
and the third term of Eq.( 179 1, due to different Kane pa- 



rameters P and different splitting A of Eg and Fy valence 
bands. Averaging these terms over the quantized motion 
in the ground state of the quantum well gives: 

8Hso - [pP^^b - A^) + 7(P2 _ Pl)] 
((|vI/.,(d/2)p-|M.,,(-d/2)|2)[kxa],- 

d 



^v^-^'-'drj_ 



^l^{^d/2)^^^,y{-d/2)) X al) (216) 



i[{nyid/2)j^^.yid/2) 
d_ 

where the constant /? is given by 

/3 



3(£;_^ 



Aw)2^ 



(217) 



scattering or due to an asymmetric potential. When only 
an asymmetric potential in the quantum well is taken into 
account, only the first term in brackets of Eq.(216l con- 
tributes to spin-orbit interactions. In order to calculate 
such term for small difference of materials parameters, 
we use the identity 

/•OO 

(V,C/,(z))=/ dz**(z)V,C/,(r)*(z)-0, (219) 



which constitutes the Ehrenfest theorem, and obtain 



|*(d/2)|2 - |*(-rf/2)|2 = --l{S/Masvm{z)). (220) 

We note that only the average of the asymmetric part 
of the confinement potential, not the full confinement 
potential, enters this expression. For gate voltage applied 
externally (S/ zUasvm(z) ) — e Ez, where E^ is the applied 
electric field. From Eq.(216), we have 



5%,, = -^[PP^{/^b - A.) +7(n' - ^^)]x 
(V^C/asy„(z))[kx a] 



(221) 



This term due to asymmetric potential defined by the 
mass and valence band splitting difference in the quan- 
tum well and barriers has the symmetry of Rashba cou- 
pling. 
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